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PREFACE 
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FORCFS  EXERTED  BY  AN  EXPANDING  ICE  SHEET 


In  trod  action 

One  of  the  most  import \nt  forces  in  the  design  of  a  dam  in  northern  regions  is  the  horizontal 
thrust  applied  to  the  structure  by  a  solid  sheet  of  ice.  Although  much  progress  has  been  made  in 
estimating  this  force,  there  is  still  a  long  way  to  go  before  we  get  an  accurate  design  value  for  it, 
mainly  because  of  the  inadequate  basic  knowledge  of  the  deformation  behavior  of  ice. 

With  time,  the  analysis  of  this  problem  has  been  more  and  more  refined  and  design  values  of 
this  pressure  have  been  successively  lowered.  At  the  beginning  of  the  century,  for  instance,  many 
dams  were  designed  for  the  condition  of  crushing  ice,1*  which  led  to  a  value  as  high  as  50  kips  per 
linear  foot. 41  At  present,  some  plastic  properties  of  ice  are  being  taken  into  account,  and  values 
used  by  Canadian  engineers  vary  between  10  and  15  kips/ft  for  rigid  structures.  For  more  flexible 
structures,  like  sluice  gates,  values  in  the  vicinity  of  5  kips/ft  are  more  commonly  used.1* 

The  importance  of  a  proper  allowance  for  ice  thrust  is  apparent  when  it  is  appreciated  that 
for  a  dam  38  ft  high,  an  ice  thrust  of  15  kips/ft  has  the  same  tendency  to  tip  the  dam  as  the  water 
pressure  it  was  designed  to  resist. 

There  have  been  few  reported  dam  failures  caused  by  ice  pressure”  ”  and  some  engineers 
have  doubted  the  advisability  of  making  substantial  allowances  for  this  pressure  in  the  design  of 
hydraulic  structures.  But  measurements  made  during  the  last  few  years  leave  no  doubt  as  to  its 
importance.  A  maximum  recorded  pressure  of  3.5  kips/ft  was  measured  on  a  tainter  gate  of  the 
Hasting  Dam14  on  the  Mississippi  River  during  the  winter  of  1932-33.  A  panel  instrumented”  and 
set  flush  with  the  face  of  Des  Joachims  Dam  during  the  winter  of  1951-52  showed  a  maximum 
pressure  of  6.9  kips/ft.  Extensive  measurements  were  made  from  1946  to  19514'  at  several 
reservoirs  in  the  mountains  of  Colorado.  Indenter  and  strain  gauges  were  set  in  thin  mortar  panels 
that  could  be  installed  on  the  face  of  the  walls  and  at  sites  exposed  to  ice  thrust.  The  maximum 
measured  ice  pressure  at  Eleven  Mile  Canyon  was  16  kips/ft  in  1947-48,  14  kips/ft  in  1948-49 
and  20  kips/ft  in  1949-50.  In  the  winter  of  1950-51  tests  were  made  at  a  number  of  reservoirs 
presenting  a  variety  of  shore  conditions.  At  Antero  and  Shadow  Mountains,  where  reservoirs  have 
flat  shores,  maximum  thrusts  of  3.6  and  5.8  kips/ft,  respectively,  were  measured.  At  the  Evergreen 
reservoir,  with  moderately  steep  shores,  the  greatest  thrust  attained  9.4  kips/ft.  The  highest 
thrust  of  all,  17  kips/ft,  was  recorded  at  Tarryall  Reservoir  which  has  steep  rocky  shores  in  the 
vicinity  of  the  dam.  Similar  tests  were  carried  out  by  VNIIG40  from  1954-56  in  the  Ges  Reservoir 
of  the  Dnieper  River.  Maximum  recorded  ice  pressure  was  8.3  kips/ft. 

The  importance  of  ice  thrust  in  the  economic  design  of  hydraulic  structures  requires  it  to  be 
evaluated  as  accurately  as  possible  at  each  site.  Its  value  depends  greatly  on  meteorological 
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conditions  and  on  tin*  restraint  caused  by  the  topography  of  the  reservoir  shores.  It  is  a  rule  of 
thumb2’  that  this  force  can  be  ignored  in  regions  where  the  mean  January  air  temperature  is  above 
10  V  (•!':■  C). 

Basic  meteorological  data 

Rising  air  temperature,  with  or  without  absorption  of  solar  radiation,  causes  an  ice  sheet  to 
expand.  The  coefficient  of  linear  expansion  of  ice'  is  close  to  28  x  10’V°F  as  shown  on  Figure  1, 
so  that  if  an  ice  sheet  1  mile  long  were  free  to  expand  it  would  lengthen  53  in.  for  a  30°F  (17°C) 
temperature  rise.  Considerable  forces  have  to  be  exerted  by  the  shores  of  a  water  body  or  by  a 
structure  to  prevent,  and  compensate  for,  this  type  of  expansion. 

It  is  not  easy  to  analyze  the  air  temperature  variations  occurring  in  nature  at  a  given  site  and 
to  choose  the  values  that  are  liable  to  produce  the  maximum  thrust  for  design  purposes.  In  most 
cold  regions  the  mean  rise  of  air  temperature  rarely  exceeds  5°F/hr  (3°C/hr)  for  a  long  enough 
interval  of  time  but  cases  have  been  known,4  mainly  in  the  Chinook  belt,*  where  the  temperature 
rose  40  F  (22"C)  in  40  minutes,  the  rise  being  followed  by  an  equally  rapid  fall  with  no  perceptible 
interval  in  between.  Many  such  variations  are  reported  in  the  literature  but  because  of  their  very 
short  duration  they  are  not  representative  of  the  temperature  fluctuation  inside  the  ice  sheet  itself 
which  lias  a  great  thermal  inertia. 

The  best  way  to  approach  the  meteorological  data  of  air  temperatures  in  a  region,  for  the 
purpose  of  obtaining  the  maximum  ice  thrust,  is  to  separate  the  temperature  rises  in  classes  with 
simple  laws  of  rise  as  a  function  of  time.  Three  types  of  simplified  analytical  laws  are  representa¬ 
tive  of  temperature  fluctu;,‘ions:  the  step  function  where  the  temperature  rises  suddenly  and  stays 
for  a  while  at  the  new  value,  the  linear  increase  for  a  long  enough  interval  of  time,  and  finally  the 
diurnal  sinusoidal  variation  of  the  air  temperature.  For  a  given  total  temperature  rise  the  first 
type  of  variation  is  the  worst  that  can  happen  to  provoke  maximum  thrust  in  an  ice  sheet.  However, 
it  cannot  be  found  at  many  places  in  nature  with  sufficiently  high  values  of  the  step  function.  The 
case  of  the  linear  rise  of  the  air  temperature  is  universal  and  can  be  found  for  extended  intervals  of 


Figure  I.  Mean  coelf  idem  of  linear  thermal  expansion  of  ice.  Standard  devia¬ 
tions  are  shown.6 


•A  region  where  strotiK.  warm,  dry  winds  of  short  duration  descend  the  eastern  slopes  of  the  Rocky  Mountains. 
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time.  For  the  same  total  air  temperature  fluctuation  the  temperature  in  the  ice  itself  is  much  more 
damped  with  a  sinusoidal  variation  than  with  a  linear  rise  and  only  this  latter  law  need  normally 
be  taken  into  account. 20 

The  analysis  of  the  characteristics  of  the  linear  rise  in  air  temperature  to  be  used  for  design 
purposes  must  take  into  account  two  variables:  the  total  temperature  rise  A  Tm  and  the  time 
interval  At  required  to  obtain  that  rise.  For  each  interval  of  time  At  from  let  us  say  5  hours  to  20 
hours  it  is  possible  to  find  the  greatest  rise  in  air  temperature  every  year  for  the  period  of  observa¬ 
tions  at  the  site.  A  typical  statistical  analysis  of  these  data20  for  23  winters  at  Quebec  City 
(Table  I)  will  then  determine  the  maximum  possible  rise  and  its  duration  for  any  expected  period  of 
return.  It  is  interesting  to  note,  for  the  case  previously  cited,  that  there  was  an  upper  limit  to 
the  rate  of  temperature  rise  (7°F/hr;  4°C/hr). 

This  case  is  not  general,  however,  and  a  predetermined  period  of  return  must  be  chosen  to 
obtain  design  values.  The  return  period  must  then  be  based  on  an  economic  valuation  criterion  to 
minimize  the  added  cost  of  construction  and  the  cost  of  the  risk  of  damage.6  This  may  essentially 
lead  to  short  periods  of  return  of  the  expected  event,  of  the  order  of  20  years. 

Table  I.  Maximum  rate  of  linear  air  temperature  rise  in  °F/hr 
for  various  time  intervals  (Quebec  City,  Canada). 


Period 


Time  intervals  in  hours 


(years) 

5 

6 

7 

8 

9 

10 

15 

20 

12.5 

4.4 

4.2 

4.3 

3.9 

3.5 

3.1 

2.1 

1.6 

25 

4.7 

4.5 

4.7 

4.1 

3.7 

3.3 

2.2 

1.7 

50 

5.2 

6.0 

5.0 

4.4 

3.9 

3.5 

2.3 

1.8 

5.6 

5.3 

5.3 

4.6 

4.1 

3.7 

2.5 

1.9 

200 

0.0 

5.7 

5.0 

4.9 

4.3 

3.9 

2.6 

2.0 

500 

6.6 

6.2 

6.8 

5.1 

4.7 

4.  1 

2.7 

2.1 

7.0 

6.5 

6. 1 

5.4 

4.8 

4.3 

2.9 

2.2 

Temperature  fluctuation  inside  the  ice  cover 

Many  factors  influence  the  temperature  distribution  inside  an  ice  cover  for  a  given  air  tempera¬ 
ture  fluctuation.  There  is  first  a  difference  between  the  temperature  at  the  ice  surface  and  the 
temperature  of  the  air  because  of  the  thermal  boundary  layer  in  the  air.  But  the  difference  is  small 
and  measurements  made  for  a  whole  winter20  on  clear  ice  have  shown  that  it  rarely  exceeds  3°F 
(1.7°C).  It  is  thus  safe  and  acceptable  to  neglect  this  effect. 

The  thermal  conductivity  of  freshly  fallen  snow  is  from  10  to  30  times  less  than  that  of  ice. 

A  few  inches  of  sncw  on  an  ice  cover  acts  as  a  practically  perfect  insulator  and  no  sensible 
temperature  fluctuation  is  transmitted  to  the  ice  sheet.  Figure  2  shows  the  dramatic  effect  of 
9  in.  of  snow  on  the  temperature  at  a  depth  of  6  in.  in  a  12-in.  ice  cover  at  Eleven  Mile  Reservoir, 
Colorado/9  This  is  one  of  the  reasons  why  young  and  relatively  thin  ice  exerts  the  highest  pres¬ 
sure  on  a  dam  at  the  beginning  of  winter.  Little  allowance  for  ice  pressure  need  be  made  during 
the  period  that  a  snow  cover  is  expected  to  stay  on  the  ice.  It  must  also  be  borne  in  mind  that 
a  uniform  fresh  snow  cover  will  also,  by  its  weight,  drown  an  ice  cover  of  approximately  the  same 
thickness,  so  that  the  mechanism  of  heat  transfer  is  completely  modified.  At  least,  no  thrust  is 
developed  while  the  ice  sheet  is  forming  again  at  the  top. 
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Figure  2.  Effect  of  a  snow  cover  on  the  temperature  at  a  depth  of  6  in.  in¬ 
side  an  ice  cover.4’ 


The  temperature  distribution  inside  an  ice  cover  can  be  obtained  from  the  Fourier  equation: 


dT 

dl 


H  * 


2  d2T 


dx‘ 


(1) 


where: 

T  -=  temperature  of  the  ice,  °F 
t  =  time,  hours 

x  =  distance  of  a  point  from  the  lower  ice  surface,  ft 
H 2  -  coefficient  of  diffusivity  of  ice,  which  can  be  taken  as  0.0434  ftVhr.1* 

The  equation  can  be  solved  for  the  simpler  boundary  conditions  discussed  previously.66  The 
simplest  case  is  that  of  determining  the  temperature  changes  (AT)  in  an  ice  sheet  originally  at  a 
steady  state  when  it  is  exposed  to  surroundings  below  freezing  temperature.  The  computation  then 
gives: 


AT 


2 


11  =  1 


An  exp 


u2  8 
Hian 

h2 


sin 


anx 


where: 


2  (sin  an  aa  cos  afl) 
%  (an  -  sin  fln  cos  an) 


and  the  an  are  the  successive  positive  roots  of  the  equation:: 


(2) 
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Here  K  is  the  thermal  conductivity  of  ice,  h  the  total  thickness  of  the  ice  sheet,  T,  the  value  of  the 
air  temperature  rise  and  Kc  the  film  conductance  of  convection  in  Btu/ftJ  hr°F. 

The  second  case  is  that  of  a  linear  temperature  rise  at  a  rate  T0.  It  then  gives: 


AT 


h2 


-  exp  - 


u2  2 

Vi. 


) 


anx 

h 


(3) 


Another  important  factor  affecting  the  rise  of  temperature  inside  an  ice  cover  is  the  absorption 
of  solar  energy.  This  effect  can  be  introduced  in  the  computation.  If  we  consider  only  the  absorbed 
heat  of  solar  radiation  at  the  surface  of  the  ice  sheet  *ji  in  Btu/ft2  hr,  we  then  obtain: 


<A  St)h 


where  Kf  is  the  film  conductance  of  radiation  in  Btu/ftJ  hr°F. 

For  the  case  of  temperature  change  due  to  absorption  within  the  ice  sheet  of  part  of  that  portion 
of  the  solar  energy  transmitted  through  the  ice,  an  additional  term  is  required  in  the  heat  equation 
which  becomes: 


d2T 
dx 2 


Tr(t)  exp  -  <A;/i 


(5) 


where  </v  is  the  absorptivity  of  ice  and  f  r  is  the  solar  volume  absorption  coefficient  of  ice  at  the 
surface  (i.e.  temperature  change  pe;  unit  time). 

The  solution  of  this  equation  for  the  proper  boundary  conditions  is: 


AT 


■{[ 


«  H'fa2  a  x' 

V'  d  i  n  ,  a 

Z  Bn  exP  -  — —  1  31tl  -J- 
n=l  h -  n 


C  ( 1  -  exp  ^‘x)  +  D~ 

r  >  h 


(6) 


where : 


2  exp  -  i p'fh 


<  h\Z 


B_ 


anil  -  cos  °n) 


a„  +  (0;*) 


[exp  \lfth  (  sin  «n  an  cos  «n)  t-  aB] 


(a„  -  sin  an  cos  a„) 


2  2 
'A; A  +  (<A;a)£ 


(Kr  +  Kc)f-  t1  "  exp  -  ^h)  (sin  -  afl  COS  aa) 


1  + 


(k,  +  I  -  -  sin  a  cos  a  \ 

\  r  °/  K  \  n  0  nJ 
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C  - 


D 


TrU)  exp  -  4r\h 

(0 p2 

T,(t)h  Tt(t)  ( Kc  +  Kf)h 

«f0;  '  Hf(0|)2 


K 


(1  -  exp  - 


1  *  (K  +  K  )  — 
v  c  x>  K 


By  means  of  the  above  equations  it  is  possible  to  compute  the  change  in  temperature  of  an  ice 
sheet  for  a  step  or  linear  function  of  temperature  rise  and  for  absorption  of  solar  radiation.  The 
heat  equation  being  linear,  the  temperature  rise  at  any  point  can  be  obtained  by  adding  up  the  effects 
of  separate  components.  This  was  done4’  for  a  40°F  (22°C)  temperature  rise  that  occurred  in  1  hour 
and  15  minutes  in  January  1947  at  Eleven  Mile  Reservoir.  The  results  are  shown  in  Figure  3.  The 
ice  temperature  is  that  of  ice  free  of  snow  at  a  depth  of  6  in.  in  an  ice  sheet  15  in.  thick.  The 
computed  values  are  shown  by  the  individual  points  near  the  smooth  observed  ice-temperature  curve. 
The  computed  and  measured  values  agree  well. 


Figure  3.  Computed  and  measured  values  of  the  ice  temperature  at  a  depth  of 

6  in.  in  the  cover.49 
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Plastic  behavior  of  ice  submitted  to  increasing  temperatures 

The  most  direct  way  of  studying  the  behavior  of  ice  under  increasing  temperature  is  to  measure 
the  force  needed  to  keep  a  test  sample  at  constant  length  between  the  loading  blocks  of  a  testing 
machine. 

The  latest  measurements  are  those  of  Monfore*'  which  were  carried  out  on  cylinders  4  in.  in 
diameter  and  4  in.  long  taken  from  lake  ice,  with  the  axes  of  the  cylinders  parallel  to  the  surface  of 
the  ice  sheet.  Loads  applied  to  the  cylinders  were  thus  in  the  same  direction,  relative  to  the  crystal 
structure,  as  thrusts  which  develop  in  an  ice  sheet  in  the  field. 

Unfortunately  Monfore  made  no  study  of  the  crystallographic  structure  of  this  ice.  Many  different 
types  of  ice  can  coexist  in  the  same  lake  and  their  creep  properties  vary  over  a  wide  range  of  values 
for  the  same  loading  conditions.  It  is  thus  now  impossible  to  ascertain  if  the  values  he  obtained  are 
representative  of  the  ice  types  that  will  give  the  highest  loading  conditions  in  nature. 

In  these  tests,  ice  temperatures  were  controlled  by  passing  air  through  the  test  chamber  and 
were  measured  inside  the  ice  with  copper-constantan  thermocouples. 

The  ice  in  the  testing  machine  was  first  brought  to  the  desired  initial  temperature  and  was  then 
held  at  this  value  until  equilibrium  was  established.  In  order  to  obtain  a  given  linear  rate  of  tempera¬ 
ture  rise  in  the  ice,  it  was  necessary  for  the  air  temperature  to  rise  at  a  much  greater  rate  for  about 
the  first  15  minutes.  As  the  temperature  of  the  ice  increased,  the  load  on  the  sample  was  adjusted 
to  maintain  the  ice  at  its  initial  length.  Temperatures  were  read  and  the  load  on  the  sample  was 
adjusted  every  5  minutes  for  the  first  30  minutes  of  test.  Temperatures  were  read  every  15  minutes 
and  the  load  adjusted  twice  every  15  minutes  thereafter. 

The  specimen  length  was  kept  constant  with  the  help  of  a  gauge  set  between  each  of  the  loading 
blocks.  The  loads  to  equilibrate  the  zero  deformation  were  applied  with  a  small  hydraulic  ram.  The 
ram  was  operated  by  a  hand  pump  equipped  with  a  large  pressure  gauge  calibrated  directly  in  total 
load.* 

More  than  100  tests  were  carried  out,  and  for  each  one  of  them  the  pressure  rose  to  a  maximum 
and  then  decreased.  Ice  temperatures  and  pressures  for  a  typical  test  are  shown  in  Figure  4.  In 
this  test  the  ice  temperature  rose  from  -10.5°F  (-23.6°C)  at  a  uniform  rate  of  5°F/hr  (2.8°C/hr). 

Check  tests  were  run  on  several  samples  to  determine  the  degree  of  reproducibility  of  apparatus 
and  techniques.  The  average  deviation  in  maximum  pressures  for  duplicate  tests  on  a  given  sample 
was  6.0%,  which  was  considered  satisfactory.  Much  larger  variations,  averaging  25%,  were  found  in 
the  maximum  pressures  reached  by  different  ice  samples  tested  under  the  same  conditions.  This 
seemB  to  depend  on  the  arrangement  of  crystals  in  apparently  similar  samples  and  it  was  eventually 
found  that  the  samples  that  reached  low  maximum  pressures  contained  crystal  surfaces  or  other 
sharp  divisions  along  which  slippage  might  have  occurred, 

Pressure  tests  were  made  with  initial  temperatures  of  -30°,  -10°  and  +10°F  (-34.4°,  -23.3°  and 
-12.2°C)  and  with  rates  of  temperature  rise  in  the  ice  of  2°,  5°,  10°  and  15°F/hr  (1.1°,  2.8°,  5.6° 
and  8.3°C/hr).  Curves  similar  to  those  of  Figure  4  were  obtained  for  each  test.  For  a  given  rate  of 
temperature  rise,  the  maximum  pressure  reached  by  the  sample  increased  as  the  initial  temperature 
decreased;  for  a  given  initial  temperature  the  maximum  pressure  increased  as  the  rate  of  temperature 
rise  increased.  The  relationships  between  the  pressures  and  time  for  various  initial  temperatures 
and  rate  of  temperature  rise  are  shown  in  Figure  5. 

All  these  results  have  basically  the  same  form  if  they  are  compared  on  a  non-dimensional 
basis  as  shown  in  Figure  6.  The  only  two  characteristic  parameters  are  the  maximum  pressure  and 
the  time  required  to  obtain  this  pressure  (Fig.  7,  8). 


•The  use  of  this  procedure  of  discontinuous  load  application  to  obtain  the  creep  behavior  of  i?e  is  very  much 
open  to  quesUon. 
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Figure  4.  Ice  temperatures  and  pressures  with  time.41 
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Figure  5.  Average  pressure-time  curves  for  a  temperature  rise  of  2°,  5n,  10°  and  15°F/hr  (1.1°, 
2.8  ,  5.6  and  8.3  C/hr)  from  indicated  initial  temperatures.41 
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Fifure  6.  Curve  of  percent  of  maximum 
pressure  and  percent  of  time  to  reach  it.4' 
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Temperature  Rise  Rota,  *f/h r 

Figure  7.  Curves  of  maximum  pressure  and 
initial  temperature  on  a  base  of  rate  of  tem¬ 
perature. rise.*1 
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a.  Creep  at  constant  stress. 


Tamperotura  Rise  Rota,  V/hr 
Figure  8.  Curves  of  time  to  reach  maximum 
pressure  and  initial  temperature  on  a  base  of 
rate  of  temperature  rise.** 


b.  Stress  at  constant  rate  of  strain. 


Figure  9.  Relationships  among  creep  parameters  for  granular-grained  ice.41 
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In  a  broad  sense  these  results  could  be  interpreted  with  the  present  knowledge  of  the  basic 
behavior  of  granular -grained  ice  under  creep  loading  conditions  for  compression  stresses  over  approx¬ 
imately  140  psi.  The  deformation  of  this  ice  as  a  function  of  time  can  be  represented41  by  curves  of 
the  type  in  Figure  9a.  In  principle  there  are  essentially  three  phases  of  ice  deformation.  At  the 
beginning  of  the  application  of  a  constant  load,  there  is  a  primary  viscoelastic  deformation.  This  is 
followed  by  secondary  deformation  at  practically  constant  rate,  the  ice  then  behaving  as  a  viscous 
substance.  In  the  final  (tertiary)  stage,  accelerated  creep  develops  up  to  the  failure  of  the  ice  sam¬ 
ple. 

It  is  also  most  interesting  to  observe  this  behavior  under  a  second  mode  of  loading.  For  a 
constant  rate  of  deformation,  the  corresponding  pressure  curve  is  of  the  form  shown  in  Figure  9b. 
During  viscoelastic  deformation,  the  pressure  rises  until  it  attains  a  maximum.  It  then  decreases 
and  stays  constant  during  the  secondary  phase  of  deformation  of  the  ice. 

There  is  little  doubt  that  this  basic  ice  behavior  is  representative  of  the  behavior  of  an  ice 
sheet  submitted  to  thermal  stresses  that  produce  significant  ice  thrusts,  as  Figure  6  shows.  The 
main  drawback  of  this  explanation  lies  in  the  fact  that  the  high  stresses  measured  by  Monfore  have 
not  been  observed  for  uniaxial  compression  with  the  values  of  strain  rates  induced  by  thermal  expan¬ 
sion.  Two  other  factors  have  also  to  be  taken  into  account  in  this  interpretation:  the  effect  of 
biaxial  expansion  and  the  influence  of  the  temperature  of  the  ice  itself  on  its  creep  properties. 

Computation  of  ice  thrust 

The  distribution  of  stresses  as  a  function  of  time  inside  an  ice  cover  subjected  to  a  linear 
temperature  rise  can  be  represented41  with  the  typical  measurements  shown  on  Figure  10.  It  can 
be  seen  that  after  the  maximum  stress  is  attained  at  the  surface  of  the  ice  sheet  the  stress  starts 
to  decrease  at  this  point  while  it  continues  to  increase  underneath.  Thus,  a  few  measurements 
show  generally  that  the  total  maximum  thrust  occurs  slightly  after  the  maximum  stress  is  attained 
at  the  upper  surface  and  its  value  may  be  from  10  to  20%  higher. 

Because  so  little  is  known  of  creep  properties  of  ice  and  because  appreciable  differences 
in  creep  rates  may  exist  between  different  kinds  of  ice  under  different  loading  conditions,  even  for 
very  small  changes  in  internal  structure,  the  analysis  was  limited  here  to  the  computation  of  the 
thrust  developed  when  the  maximum  stress  was  attained  at  the  upper  surface  for  Monfore’s  ice.  At 
the  present  time  this  is  all  that  can  be  used  as  an  indication  for  design  purposes  and  it  is  expected 
to  give  values  very  much  on  the  high  side. 

The  total  thrust  was  computed  with  the  usual  assumptions  for  this  problem:40 

1.  The  initial  air  temperature  is  -30°F,  -10° F  and  10°F  (-34.4°,  -23. 3n  and  -12.2°C). 

2.  The  initial  ice  temperature  varies  linearly  from  -30°F  (-34.4°C)  at  the  upper  surface 
to  *32  'F  (0  C)  at  the  lower  surface. 

3.  The  air  temperature  changes  at  rates  of  2°  to  15°F/hr  (1.1°  to  8.3°C/hr). 

4.  The  temperature  of  the  ice  at  the  surface  is  assumed  to  oe  the  same  as  the  air 
temperature. 

5.  The  thickness  of  the  ice  sheet  remains  constant. 

The  temperature  distribution  and  the  rate  of  temperature  rise  inside  the  ice  sheet  were  computed” 
by  eq  3  for  variable  ice  sheet  thicknesses,  air  temperature  rises  and  initial  ice  temperatures. 

Figure  11  shows  that  for  an  ice  sheet  of  given  thickness  the  maximum  thrust  is  attained  for  an 
air  temperature  rise  of  only  3°F/hr  (1.7°C/hr).  This  can  be  explained  by  the  fact  that  at  this  rate 
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PRESSURE.  p«< 


Figure  10.  Measured  pressure  distribution  in  a  cross  section  of  an  ice  sheet." 
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Without  solar  energy  and  restraint 


Pate  Of  Rise  Of  Air  T«mp*rotur»,*F/hr 


Figure  11.  Thrust  of  an  ice  sheet  and  rate  of  rise  of  air 
temperature. J0 


there  is  enough  time  for  the  ice  to  develop  an  important  overall  thrust  before  the  maximum  pressure 
is  attained  at  the  surface  when  the  computation  is  stopped.  The  curves  also  show  an  inflection 
point  at  a  rate  of  about  7.6°F/hr  (4.2°C/hr).  This  can  be  explained  from  the  rheological  behavior 
of  Monfore's  ice  in  this  area  (Fig.  7)  where  there  is  a  sharp  change  in  the  variation  of  the  maximum 
pressure  that  this  ice  can  exert. 
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Ice  Sheet  Thickness,  inches 

Figure  12.  Maximum  thrust  exerted  by  an  expanding 
ice  sheet  lor  various  ice  thicknesses  at  three  latitudes 
and  three  initial  air  temperatures. " 


r  effect  of  ice  thickness  is  shown  in  Figure  12.  The  thrust  does  not  increase  appreciably 
or  tee  sheets  thicker  than  12  in.,  mainly  because  there  is  no  sensible  variation  of  temperature 
*\ow  this  point  before  the  stresses  start  to  decrease  at  the  upper  surface.  Furthermore  if  we 
also  consider  that  snow  covers  are  related  to  thicker  ice  sheets,  there  is  little  doubt  that  the 
maximum  thrust  does  not  increase  appreciably  for  ice  sheets  thicker  than  12  in. 

When  the  ice  sheet  is  free  of  snow  and  exposed  to  sunlight,  its  temperature  may  have  an  added 

increment  by  absorption  of  solar  radiation,  quite  apart  from  the  change  caused  by  an  increase  in 
air  temperature. 

The  magnitudes  of  ice  thrusts  were  computed  to  include  solar  energy  effects  with  the  following 
additional  assumpMons: K  B 

The  latitudes  are  40°,  50°  and  60°. 

The  vernal  equinox  is  the  time  limit  for  the  formation  of  clear  solid  ice  of  adequate 


t . 
8. 

thickness. 


9.  The  atmospheric  transmission 
high  elevation. 


constant  is  0.9,  based  on  a  clear  atmosphere  at  a  fairly 


10.  The  coefficient  of  heat  transfer  is  2  Btu/ftJ  hr  °F, 
absorbed  solar  energy. 


to  allow  for  surface  losses  of  the 


For  these  assumed  conditions  the  solar  energy  that  reaches  the  surface  of  the  ice  sheet  durinc 
nudday  is  about  250  Btu/ft1  hr  at  40 1  latitude.  200  Btu/ftJ  hr  at  50°  latitude  and  150  Btu/ftJ  hr  at 
60  latitude  The  absorptivity  of  ice  is  small  in  the  visible  spectrum”  but  beyond  a  wavelength  of 
micron  it  is  great  and  is  analogous  to  that  of  water.  Consequently,  about  40%  of  the  solar  energy 
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reaching  the  surface,  essentially  all  taat  of  wavelength  greater  than  1  micron,  is  absorbed  within 
4  in.  of  the  surface.  For  a  clear  ice  sheet  1  ft  thick  an  additional  15%  of  the  energy  is  absorbed, 
whereas  for  a  4-ft  thickness  this  ad  Jitional  absorption  amounts  to  about  25%.  Most  of  the  remainder 
of  the  energy  is  transmitted  although  a  small  percentage  is  reflected. 

The  effect  of  solar  energy  was  computed"  ucing  eq  4  find  the  results  (Fig.  12)  indicate  that  the 
absorption  of  this  energy  may  produce  a  considerable  increase  in  thrust  above  that  caused  by  the 
rise  of  the  air  temperature. 

The  effect  of  lateral  restraint  was  estimated  on  the  basis  of  the  behavior  of  an  elastic  slab. 
Elastic  behavior  is  approached  only  for  the  conditions  of  very  rapid  strain  rates  when  flow  is  small. 
Poisson’s  ratio  v  may  not  be  applicable  for  flow  conditions,  but  pressures  based  on  elastic  be¬ 
havior  for  restraint  conditions  might  give  the  most  severe  loading  conditions.  The  ratio  of  inciease 
in  stress  for  complete  lateral  restraint  to  the  stress  developed  for  no  lateral  restraint  is  given  by 
iV(l  -  u).  The  values  computed  for  no  lateral  restraint  are  increased  by  this  factor  and  are  shown 
in  Table  II. 


Table  n.  Total  thrust  exerted  by  expanding  ice  sheet  (Monfore's  ice). 
Ice  thickness  greater  than  12  in.,  no  solar  energy. 


Initial  air  temperature 

S0°F 

-10°  F 

+10°  F 

Ice  aurtace 

Lateral  restraint 

Lateral  restraint 

Lateral  restraint 

temp  rtae 

No 

Fes 

No 

Yes 

No 

Fes 

S°F/hr 

11.6 

18 

a.t 

ie.7 

4.6 

7.1 

6°F/hr 

10.0 

16.6 

7.3 

11.3 

3.5 

5.5 

10°  F/hi 

8.8 

13.8 

6.8 

9.1 

2.4 

3.7 

16°F/hr 

8.6 

13.4 

5.0 

8.9 

2.2 

3.5 

Secondary  factors 

Observations  on  lakes  free  of  snow40  have  shown  that  reefs  and  pressure  ridges  formed  by  the 
expansion  and  buckling  of  an  ice  sheet  seldom  occur  when  the  ice  is  thicker  than  1  ft,  although 
they  have  been  known  to  form  in  ice  20  in.  thick.  This  phenomenon  might  be  interpreted  to  indicate 
that  for  ice  sheets  less  than  1  ft  thick  the  critical  buckling  load  is  the  limiting  factor  governing 
the  ice  pressure  on  a  dam,  whereas  for  thick  ice  sheets  the  pressure  may  develop  to  its  maximum 
without  buckling  of  the  sheet.  These  observations  have  been  somewhat  correlated  by  a  computa¬ 
tion  in  the  first  approximation,41  assuming  very  approximate  plastic  ice  properties,  where  it  is  shown 
that  buckling  is  a  limiting  factor  when  the  ice  thickness  is  less  than  1.5  ft  thick.  For  ice  sheets 
2  ft  thick,  which  would  normally  be  considered  for  design  purposes,  this  limitation  need  not  be  taken 
into  account. 

When  an  ice  sheet  increases  in  thickness  in  a  reservoir  and  the  temperature  falls,  tensile 
stresses  are  Bet  up  in  the  ice,  giving  rise  to  contraction  cracks.  When  the  ice  is  heated  up  again 
the  gap  in  the  crack,  if  it  has  not  already  frozen,  will  first  have  to  be  taken  up  by  the  expanding 
ice  before  any  thrust  can  be  exerted.  This  phenomenon  thus  reduces  the  force  exerted  by  an 
expanding  ice  sheet. 

Many  other  factors  tend  to  reduce  the  thrust  of  ice  compared  to  the  thrust  that  can  be  computed. 
An  important  one  is  the  difference  between  the  behavior  of  a  small  ice  sample  in  a  laboratory  and 


M 


ICE  PRESSURE  ON  ENGINEERING  STRUCTURES 


a  large  ice  sheet  in  nature.  Because  of  cracks,  faults,  internal  dislocations  and  discontinuity  in  a 
lar^e  ice  sheet  a  small  sample  is  always  stronger  than  the  ice  in  the  field.  It  has  also  been  shown11 
that  the  creep  rate  of  ice  increases  under  reloading  conditions.  An  ice  sheet  submitted  to  cycles  of 
contraction  and  expansion  during  its  growth  will  exert  much  lower  pressure  than  what  would  be 
expected  fiom  the  belt.  '  lor  of  a  laboratory  sample  tested  forthe  first  time. 

Static  pressure  against  Isolated  structures 

The  pressure  of  an  ex  >anding  ice  sheet  on  the  tip  of  a  protruding  spillway  pier  set  in  the  ice  is 
somewhat  higher  than  the  pressure  per  unit  length  multiplied  bv  the  width  of  the  pier.  A  first  approxi¬ 
mation  of  this  value  eai  '  i  obtained  with  an  empirical  formula:’® 

f  p„„  (»o  •  ~)  (7) 

where: 

F  maximum  horizontal  force,  lb 

p  static  pressure  per  unit  length  obtained  for  this  case,  lb/ft 

width  of  the  pier,  ft 
B'  span  between  piers,  ft. 

For  the  case  of  a  bridge  pier  this  force  has  to  be  taken  into  account  only  if  it  is  not  balanced 
by  an  equal  and  opposite  force  on  the  downstream  side  of  the  pier.  Because,  in  a  general  way, 
the  turbulence  of  the  flow  in  the  wake  of  a  pier  might  prevent  ice  from  forming  close  to  it,  it 
might  be  safer  in  most  applications  to  consider  this  force  on  the  upstream  edge. 

Furthermore,  the  lateral  ice  pressures  on  a  bridge  pier  usually  balance  each  other  on  both 
sides.  There  are  many  cases,  however,  in  high-velocity  flow,  where  a  channel  may  stay  open  in  the 
winter  between  two  bridge  piers.  If  an  ice  sheet  then  forms  on  the  other  side  of  these  particular 
piers,  full  allowance  has  to  be  made  for  lateral  pressure.  Such  a  state  of  loading  may  greatly  affect 
the  stability  of  piers. 

Protection  against  static  ice  pressure 

At  an  existing  dam,  designed  with  no  allowance  for  ice  pressure,  the  following  procedure  may 
be  used  to  protect  against  expanding  ice.  Generally,  all  dams  are  designed  for  the  full  reservoir 
condition  with  the  water  surface  at  or  very  near  the  top  of  the  dam.  It  is  possible  to  compensate 
for  the  effect  ol  ice  thrust  by  lowering  the  reservoir  level  during  the  season  when  ice  pressure  may 
exist,  so  that  the  total  load  at  the  critical  section  is  not  larger  than  for  the  full  reservoir  condition. 

Figure  13®°  shows  the  required  lowering  of  a  reservoir  to  balance  the  effect  of  the  indicated 
ice  thrusts.  The  criterion  used  for  evaluation  in  Figure  13a  is  the  moment  at  the  base,  whereas 
Figure  13b  shows  results  with  horizontal  shear  as  the  governing  factor.  The  moment  effeci  governs 
(that  is,  yields  larger  reductions  in  the  water  level)  for  all  heigh® s  greater  than  those  at  the  maxi¬ 
mum  points  of  the  curves  in  Figure  13a.  The  height  above  which  the  moment  governs  and  below 
which  shear  governs  is  shown  by  the  lines  marked  "moment-shear  line."  The  base  of  a  dam  is 
usually  the  critical  section,  but  the  curves  can  be  applied  to  any  other  elevation  by  considering  the 
height  above  the  elevation  in  question.  The  curves  are  applicable  to  a  gravity-type  dam. 

The  curves  indicate  that  in  the  case  of  high  dams  the  ice  thrust  is  of  relatively  minor  impor¬ 
tance  in  its  effects  on  forces  at  the  base  of  the  section,  whereas  for  low  dams  it  is  of  considerable 
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Figure  13.  Reduction  in  reservoir  water  level  to  compensate  lor  ice  thrust. 60 

importance.  For  high  dams,  it  should  be  borne  in  mind  that  an  elevation  at,  say,  mid-height  might 
give  a  condition  more  severe  than  at  the  base  elevation,  even  though  the  latter  is  the  location  ol 
maximum  stress  for  full  reservoir.  Complete  consideration  of  this  factor  requires  actual  stress  and 
stability  analyses. 

The  best-known  method  of  preventing  an  ice  sheet  from  pushing  on  a  dam  is  to  maintain  an 
open  water  gap  in  front  of  the  stiucture.  This  can  be  done  in  various  ways.  A  permanent  solution 
is  to  install  a  system  of  piping  along  the  edge  of  the  wall  at  the  bottom  of  the  reservoir  so  that 
compressed  air  can  be  used  to  create  a  vertical  current  bringing  the  warm  water  at  the  reservoir 
bottom  to  the  top,  thus  preventing  ice  formation.  This  system  works  well  as  long  as  there  is  an 
adequate  vertical  temperature  gradient  in  the  reservoir.  Ir  turbulent  flows  this  gradient  is  very 
small  and  there  is  normally  not  enough  heat  stored  in  the  bottom  layer  of  water  to  prevent  ice 
formation  at  the  top. 

A  wooden  boom  laid  along  the  perimeter  of  the  wall  of  a  reservoir  has  been  found  very  effective 
in  protecting  it  againBt  ice  pressure.17  Although  no  complete  explanation  was  ever  given  for  this 
beneficial  effect  it  seems  that  the  boom  actB  as  a  hinge  when  the  ice  rises  and  falls  and  as  a  roller 
when  it  expands  and  contracts.  Another  solution,  which  has  never  been  tried,  would  be  the  setting 
of  plastic  foam,  or  rubber,  panels  at  the  water  line  on  the  walls  of  the  structure.  If  thick  enough, 
these  materials  could  take  a  good  amount  of  deformation  without  transmitting  any  sensible  thrust 
to  the  wall. 

In  an  emergency  a  gap  a  few  inches  wide  can  be  made  by  hand  with  ice  chisels,  saws  or  steam 
heating.  Chemicals  might  also  be  used  to  melt  and  considerably  weaken  the  ice  in  front  of  the 
endangered  structure.  No  unassisted  relief  to  the  situation  can  be  expected  before  an  adequate  snow 
cover  has  formed  on  the  ice  sheet. 
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Shore  effects  la  thermal  expansion  of  lake  ice 

Expansion  and  contraction  of  lake  ice  due  to  temperature  changes  has  long  been  known  to  be 
responsible  for  the  formation  of  ice  ramparts  along  the  shore  areas,  or  of  pressure  ridges  In  the  ice 
itself. 

The  underlying  principle  involved  in  the  formation  of  ice  ramparts  has  been  well  established.1 
A  rapid  fall  in  air  temperature  causes  the  ice  on  a  lake  to  contract,  so  that  tension  cracks  are 
developed  in  the  ice  sheet.  These  cracks  immediately  fill  with  water,  which  quickly  flreeaes,  so  that 
the  total  mass  of  the  ice  cover  is  increased.  A  subsequent  rise  in  air  temperature  results  in  expansion 
of  the  ice,  so  that  the  total  surface  area  is  greater  than  it  was  before  contraction.  This  increase  in 
area  of  the  ice  cover  results  in  compressive  forces  against  the  shore  so  that  shore  debris  is  forced  into 
ridges  called  "ice  ramparts"  (Fig.  14). 

Extensive  measurements  were  made  of  this  phenomenon  at  Wamplers  Lake,  Michigan,  during  ihe 
winters  of  1951-52  and  1952-53.”  The  first  winter  a  movement  of  ice  on  the  shore  of  1.6  ft  was 
measured,  of  the  type  shown  in  Figure  15.  The  expansion  movement  of  the  ice  was  not  normal  to 
the  shore  but  was  related  to  the  form  of  the  lake;  more  expansion  was  observed  along  the  long  axis 
than  along  the  short  one  on  this  oblique  lake.  In  the  second  winter  the  lake  level  was  lower  and  a 
new  rampart  was  built  at  a  lower  level  from  boulders  and  gravel  that  were  normally  submerged  be¬ 
neath  the  water  of  the  shore  region. 

On  one  occasion  the  atmospheric  temperature  dropped  from  30°  to  8°F  (-1°  to  -13°C)  in  12  hours 
and  cracks  developed  with  widths  ranging  from  a  few  tenths  of  an  inch  to  1.5  in.  The  ice  thickness 
at  that  time  was  8  in.  and  the  snow  cover  was  3 Vi  in.  One  set  of  cracks  radiated  from  the  central 
part  of  the  lake  and  another  set  roughly  followed  the  shoreline.  Figure  46  shows  two  of  these 
orthogonal  cracks. 

It  is  interesting  that  the  formation  of  a  compression  ridge  was  observed  only  once  during  these 
two  winters,  when  the  ice  was  8  in.  thick.  The  axis  of  this  ridge  lay  at  right  angles  to  the  direction 
of  movement  of  the  sheet  as  recorded  by  markers.  Prior  to  this  time  no  buckling  had  taken  place, 
but  considerable  thrust  had  developed  on  the  shore. 


Figure  14.  Ice  ramparts.” 
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Figure  16.  Orthogonal  cracks  on  Wamplers  Lake,  Michigan  71 
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Finally  it  was  observed  that  no  appreciable  movement  of  the  markers  occurred  because  of  the 
skin  friction  caused  by  the  wind. 


IMPACT  FORCES  OF  ICE  ON  8TRUCTURE8 


IaUodactlon 

Most  older  river  bridges  of  northern  countries  end  their  lives  by  pier  movement  following  scour 
and  ice  pressure  at  breakup  time.  Very  little  has  been  written  on  this  unfortunate  event,  although 
it  would  be  quite  revealing  for  the  design  of  these  structures.  There  is  no  doubt  that  one  of  the 
most  important  forces  that  can  lead  to  the  collapse  of  a  bridge  is  the  force  of  moving  ice  floes  but 
only  the  failure  of  larger  bridges  under  ice  conditions  has  been,  at  times,  reported  in  the  literature. 
Examples  are  the  Niagara  Falls  View  bridge,"  the  Cushing  bridge,"  some  Missouri  highway  bridges," 
and  the  Brattleboro  bridge.11 

Ice  may  adversely  affect  a  bridge  in  a  number  of  ways.  Premature  snow  melt  in  a  river  basin 
may  suddenly  increase  the  river  flow,  which  in  turn  raises  the  solid  ice  sheet  adhering  to  the  piers, 
partially  lifting  them  out  and  endangering  the  stability  of  their  foundation.  The  piers  may  be  located 
in  such  a  manner  in  the  river  bed  as  to  hinder  the  passage  of  ice,  and  to  prime  ice  jams.  These 
increase  the  pressure  on  the  piers,  and  the  superstructure  may  be  destroyed  if  the  ice  is  raised  high 
enough  to  impinge  on  it.  The  water  may  also  find  a  new  channel  within  or  around  the  ice  jam  and 
high-velocity  flow  may  badly  scolir  the  approaches  or  the  pier  foundation.  Finally  the  moving  ice 
floes  at  breakup  usually  exert  an  appreciable  force  when  hitting  the  bridge  piers.  This  force  may 
displace  a  pier  enough  to  cause  the  collapse  of  the  bridge. 

Moving  ice  may  also  impinge  on  other  types  of  hydraulic  structures.  Figure  17  shows  a  steel 
dam  gate  that  failed  under  the  impact  of  ice. 


Figure  17.  Steel  stop-log,  88  ft  long,  broken  by  ice  impact  (1967). 
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ui  North  America  the  practice  has  generally  been  to  compute  the  impact  force  of  ice  on  struc¬ 
tures  a  1th  the  design  code  of  the  AASHO.1  It  simply  says  that  one  has  to  apply  a  horizontal  com¬ 
ponent  of  force  to  the  pier  equal  to  the  projection  of  the  ice  area  perpendicular  to  the  axis  of  the 
flow,  multiplied  by  the  crushing  strength  of  ice,  which  is  recommended  to  be  taken  as  400  psi.  This 
amputation  usually  leads  to  considerable  ice  forces  that  increase  the  cost  of  bridges.  The  code 
values  are  widely  questioned  by  bridge  designers  because  of  the  supposedly  high  value  taken  for 
the  crushing  strength  of  ice  and  mainly  because  they  do  not  take  into  account  any  mode  of  failure 
Other  thau  compression  although  piers  are  designed  specially  to  shear  or  lift  the  floes  out  of  the 
water  in  order  to  reduce  this  force.  We  will  now  review  the  most  important  recent  work  on  this  ques¬ 
tion,  principally  that  of  Korzhavin.40 

Strength  of  Ice  nnder  impact  loading 

Under  conditions  of  quick  loading  as  when  an  ice  floe  hits  a  bridge  pier,  ice  behaves  like  a 
brittle  material  and  its  ultimate  strength  can  be  determined.  Because  ice  is  formed  in  nature  under 
various  conditions  there  is  a  wide  spread  in  the  results  of  tests  made  to  determine  its  mechanical 
properties.  This  question  is  discussed  extensively  in  the  CRREL  monographs  on  physics  and 
mechanics  of  ice  (see  CRSE  Monograph  D-C2,  Ice  as  a  material *)  and  we  will  discuss  here  the 
engineering  properties  of  interest  for  this  application. 

We  are  concerned  with  river  and  lake  ice  which  can  usually  be  divided  into  clear  lake  ice  and 
snow  ice.  The  clear  lake  ice  usually  consists  either  of  a  fabric  of  prisms  with  their  lengths  normal 
to  the  ice  sheet,  or  of  oongealed  frazil  slush.  The  columnal  crystals  usually  increase  in  diameter 
with  depth  and  the  number  of  crystals  at  a  given  depth  decreases  accordingly,  rhe  crystal  size 
varies,  depending  probably  on  the  freezing  conditions,  from  a  fraction  of  an  inch  to  several  inches 
across.  On  the  surface  of  the  sheet  there  is  usually  a  layer  of  ice,  called  snow-ice,  made  of  small 
compact  grains.  In  a  lake  this  ice  is  formed  either  from  a  snowfall  freezing  into  the  water  surface 
or  when  water  sweeps  through  the  ice  cover,  saturates  the  snow  cover,  and  then  freezes.  The 
freezing  of  slush  and  frazil  forms  congealed  frazil  slush  which  is  normally  transparent  ice  of 
granular  structure.  Snow  ice  usually  contains  many  air  bubbles  and  is  weaker  than  'jlear  lake  ice. 

The  mechanical  properties  of  river  and  lake  ice  depend  on  a  number  of  factors,  the  most 
important  being  the  temperature  of  the  ice,  its  internal  structure  and  the  conditions  of  loading.  Be¬ 
cause  of  these  factors  there  is  a  wide  variation  in  ice  strength  measurements  quoted  in  the  litera¬ 
ture  for  thousands  of  tests  on  lake  ice.  The  extreme  range  of  results  is  shown  in  Table  ID. 

Table  m.  Range  of  Ice  strength  as  Massed. 

Compressive  strength  60  to  1800  psi 

Flexural  strength  44  to  414  psi 

Shear  strength  49  to  182  pel 


The  temperature  of  fee  is  one  of  the  most  important  factors  affecting  its  strength.  Figure  18 
gives  the  crushing  strengths  of  strong  clear  lake  ice  and  snow  ice.7  The  lake  ice  was  columnar, 
made  of  prismatic  crystals  to  2  cm  by  4  to  7  cm  with  the  crystallographic  c-axis  nearly  parallel  to 
the  surface  of  the  sheet.  The  snow  ice  bad  Bmall  granular  crystals  tt  to  2  cm  in  diameter  with  many 
air  bubbles  and  a  density  from  0.82  to  0.87  g/cm\  Ten  tests  were  made  at  each  temperature.  It  can 
be  seen  that  for  a  load  applied  parallel  to  the  ice  sheet  with  clear  lake  ice  the  strength  was  177  psi 
at  32°  F  (0°C),  400  psi  at  23°  F  (-6°C)  and  went  as  high  as  1000  psi  at  -58°F  (-50°C). 


•In  preparation. 
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a.  Clear  lake  ice. 
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Figure  18.  Crushing  strength  vs  temperature  of  ice  on  Portage  Lake,  Michigan.7 


Flit*  structure  of  ice  also  has  a  great  influence  on  its  strength.  The  quality  of  the  melt  is 
important  because  salts,  gases  and  organic  impurities  sometimes  collect,  with  air,  at  the  boundaries 
of  grains  during  their  formation.  The  crystal  structure  -  its  form,  dimensions  and  arrangement  and 
the  relative  axis  orientation  of  the  grains  -  characterizes  types  of  ice  of  widely  different  properties. 
The  general  faults  in  the  structure  such  as  stratification,  cracks  due  to  previous  loading  and  thermal 
conditions,  age  and  grain  boundary  melt  (candling)  play  a  dominant  role  in  ice  strength.  This  last 
factor  bears  particularly  on  the  present  application.  It  was  found3'  that  a  combination  of  bright  sun 
and  air  and  ice  temperatures  near  32  F  (0°C)  will  decrease  the  bending  strength  of  ice  by  as  much 
as  six  times  from  morning  to  midafternoon. 

The  strength  of  ice  also  depends  on  the  conditions  under  which  it  is  tested.  The  direction  of 
stress  application  relative  to  crystal  structure  and  the  size  of  the  specimens  have  an  influence.9 
The  flexural  strength  of  small  beams  is  close  to  twice  as  high  as  that  of  bigger  beams  at  the  same 
location.37  The  rate  of  application  and  the  duration  of  the  lead  are  also  very  important.  Figure  19 
shows  the  results  of  tests  to  determine  the  influence  of  the  strain  rate  on  the  crushing  strength  for 
columnar  ice  with  horizontal  c-axes.13 

To  decide  upon  a  value  ior  the  strength  of  ice  in  the  design  of  isolated  structures  in  lakes  and 
rivers  two  different  conditions  may  have  to  he  considered.  The  first  is  the  movement  of  ice  floes  by 
current  and  wind  during  the  wilder,  especially  m  lakes  where  the  ice  sheet  does  not  form  over  the 
whole  surface.  Under  these  conditions  the  temperature  of  the  ice  is,  on  the  average,  a  few  degiees 
tie  low  32  F  i0  C).  It  must  he  32  at  the  ice  water  interface  and  even  in  very  cold  weather  it  does 
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Figure  19.  Influence  ol  strain  rate  on  the  crushing  strength  of  laboratory- 
made  columnar  ice  with  horizontal  c-axes.  Temperature :  -10° F  (-23°C).12 


not  usually  get  much  below  32°  at  the  surface  because  of  the  insulating  effect  of  the  snow  cover. 
Besides,  the  temperature  distribution  is  far  from  linear  in  the  ice  sheet  itself,  and  the  tempera¬ 
ture  in  the  upper  layer  is  quickly  attenuated  toward  32°  at  the  lower  interface. 

The  second  condition,  the  impact  of  ice  floes  at  spring  breakup,  is  the  usual  condition  for  the 
design  of  a  bridge  pier  in  a  river.  The  water  has  previously  flowed  over  the  ice  sheet,  which  has 
been  considerably  weakened  by  solar  radiation  before  breaking  up  and  moving  down  the  river. 

For  these  two  conditions  of  winter  and  spring  ice  the  strength  values  are  chosen  according  to 
existing  design  values  as  given  in  Table  IV  for  example. 

Table  IV.  Design  values  of  ice  strength  for  impact  loading. 


Winter  ice 

Spring  ice 

Compressive  strength  <Jq 

400  psi 

75  psi 

Flexural  strength  <7^ 

200  psi 

65  psi 

Shear  strength  rQ 

120  psi 

60  psi 

For  winter  ice  strength  design  values  should  be  taken  for  the  strongest  type  of  ice  that  can  be 
produced  in  sizable  area  in  a  laky  or  a  river.  Unfortunately  this  is  unknown  at  the  present  time  but 
for  clear  lake  ice  with  big  prismatic  crystals  and  little  occluded  foreign  substances  values  of  aQ  = 
400  psi  and  rQ  =  120  psi  have  been  found  at  23° F  (-5°C).7  This  crushing  strength  value  of  400  psi 
is  the  same  as  that  suggested  by  a  Committee  of  the  ASCE  in  1931“  and  introduced  in  the  specifica¬ 
tions  of  the  AASHO,  The  strength  of  ice  in  bending,  Table  IV,  is  suggested  for  small  specimens  of 
clear  lake  ice  at  23°F.M  Before  more  is  known  about  the  strength  of  river  and  lake  ice  this  may  be 
used  for  design  practice  for  a  solid  clear  columnar  ice  sheet  where  by  chance  configuration  no 
sensible  flaw  is  present  in  the  area  under  impact  loading.  Actual  testing  in  nature  with  big  ice 
beams  shows  a  flexural  strength  which  is  normally  lower.’  "  21 
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The  strength  values  for  spring  ice  were  taken  directly  from  the  highest  values  recommended  by 
Korzhavin  for  rivers  of  the  Northern  U.S.S.R.  and  Siberia.40  This  is  for  ice  floes  moving  at  a  speed 
of  3  ft/sec  and  the  diiection  of  the  force  essentially  normal  to  the  crystal  prismatic  axis.  These 
values  were  obtained  after  standard  testing  and  also  by  measurements  of  forces  acting  on  piers. 

Indentation  of  Ice  by  a  pier  with  vertical  udge 

The  destruction  of  an  ice  floe  by  a  pier  with  a  vertical  edge  occurs  in  different  ways,  depending 
mostly  on  the  stored  kinetic  energy  of  the  floe  and  its  strength.  Upon  impact  with  a  pier  a  small  ice 
floe  is  locally  crushed,  stops  and  then  continues  to  move  across  an  adjacent  opening.  A  slightly 
larger  floe  generally  collapses  by  splitting  while  a  large  ice  field  is  cut  without  splitting.  During 
the  penetration  of  a  pier  into  a  large  ice  field,  the  width  of  indentation  gradually  increases  and 
attains  a  maximum  value  equal  to  the  total  widtli  of  the  pier,  after  which  it  remains  constant. 

General  relationships.  If  we  consider  an  ice  floe  hitting  a  pier,  the  energy  transmitted  to  the 
obstacle  is  made  up  of  the  energy  of  deformation  inside  the  ice  floe  and  the  energy  required  for  its 
crushing.  As  the  energy  of  deformation  is  negligible  compared  with  the  energy  of  crushing  for 
sizable  impacts,  the  fundamental  relationship  of  mechanics  gives: 

F  dx  m  v  dv  (8) 

where  F  is  the  force  exerted  on  the  obstacle  for  a  displacement  dx  of  the  floe,  and  m  and  v  are  the 
mass  and  instantaneous  velocity  of  the  floe. 

On  the  other  side,  at  crushing,  we  have: 

F  o(  y  h  (9) 

where  <»  is  the  crushing  strength  of  ice  by  indentation,  y  is  the  width  of  the  failure  zone  and  h  is 
the  thickness  of  the  ice. 

For  an  ice  floe,  with  an  initial  velocity  F0,  that  is  completely  stopped  by  the  pier  we  get,  using 
eq  8  and  9: 

x  xQ 

f  h  ui  (x)  y  (x)  dx  (m  Vq).  (10) 

x  o  2 


Because  of  the  geometry  of  the  pier,  the  term  on  the  left  has  an  upper  limit  when  x  dQ  and 
y  B0  where  dQ  is  the  length  of  the  pier  nose  and  B0  is  the  total  width  of  the  pier. 

If  the  kinetic  energy  in  eq  10  is  adequate  to  attain  this  upper  limit,  the  floe  will  not  stop  and 
the  maximum  force  will  be  given  by: 


",  B0h. 


(ID 


In  the  case  where  the  indentation  strength  of  the  ice  can  be  taken  as  independent  of  the  indenter 
width  x,  the  force  required  to  stop  an  ice  floe  of  a  limited  size  can  be  easily  computed: 
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yo 

o 


y  (x)  dx 
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Figure  20.  Crushing  ol  an  ice  lloe  on  a  pier. 


(12) 


H  yQ  h  . 


Aq  is  thus  the  horizontal  crushed  surface  of 
the  ice  floe  when  it  stops  as  shown  on  Figure  20 
and  can  be  determined  for  any  type  of  pier  shape." 
Actually,  y0  is  the  pier  width  corresponding  to  the 
Aq  area. 


For  the  case  of  a  triangular-nosed  pier  of  angle  2a  attacked  by  an  ice  floe  along  its  axis  we 
obtain  the  Pctrunichev  formula"  for  limited  size  floes: 


Hmax  =  V0  h  \/ 2  P° i  atan  a  (13) 

where: 

=  maximum  horizontal  force,  lb 

111  d  A 

F0  =  water  velocity,  ft/sec 
h  =  ice  thickness,  ft 
p  =  density  of  ice,  1.78  lbf-sec2/ft4 
ft  =  area  of  the  ice  floe,  ft2 
cr  =  indentation  strength  of  ice,  lb/ft2  . 

Strength  ol  ice  against  ind?ntation.  The  effect  of  the  indenter  dimensions  on  the  strength  of 
ice  was  studied  experimentally  by  Korzhavin/0  The  ice  used  in  the  tests  was  taken  from  the  river 
Ob  at  Novosibirsk.  Its  compressive  strength  was  determined  on  22  samples  70  x  70  x  70  mm  at 
32° F  (0°C).  The  effect  of  load  indentation  was  studied  on  177  samples  with  sizes  varying  between 
70  x  70  x  200  mm  and  70  x  70  x  300  mm.  The  tests  were  performed  as  shown  in  Figure  21  and 
the  following  empirical  relationship  was  proposed: 

■  °°  VsT  (14) 


where  oQ  and  a t  are  respectively  the  strength  of  ice  in  uniaxial  compression  and  its  strength  with 
an  indenter  of  width  bQ  acting  on  an  ice  sample  of  length  B.  There  is  a  limit  to  the  value  of  B 
that  can  be  used  in  the  formula.  For  a  large  ice  floe  a  complicated  state  of  stress  appears. 
Immediately  following  contact  with  the  pier  there  is  a  zone  of  crushed  material  surrounded  by  zones 
of  plastic,  and  elastic,  deformation  as  shown  on  Figure  22.  Taking  the  stress  distribution  of 
Boussinesq  in  the  elastic  zone,  Korzhavin  considers  that  the  practical  limit  of  induced  stresses  in 
the  ice  floe  is  given  by  the  isobar  corresponding  to  about  2  to  3%  of  the  mean  indenting  pressure. 
The  width  B{  of  this  isobar  is  approximately  15  ft0,  which  can  be  taken  as  the  effective  width  of 
large  ice  floes  (B  >  B  j).  Consequently  when  a  large  ice  floe  is  cut  by  a  pier  the  indentation 
strength  cannot  exceed  the  value: 
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Figure  21.  Indentation  test. 


Figure  22.  Stress  zones  in  an  indented  block  of  ice. 


(15) 


Influence  of  the  rate  of  strain  on  the  strength  of  ice.  Tests  carried  out  by  G.A.  Pchelkin, 

N.A.  Tsytovtch,  and  S.S.  Vialov  show  that  the  strength  of  ice  decreases  with  increasing  strain  rates. 
Korzhavin  and  I’utkhin40  made  a  systematic  study  of  this  factor.  The  range  of  deformation  studied 
on  cube  samples  with  10-cm  sides  was  from  2  to  32  cm/min.  They  found  that  the  variation  of 
strength  of  ice  for  rates  of  strain  from  10'1  to  10'Vsec  could  be  expressed  by  the  following  empirical 
relation 


where  I  is  the  axial  strain  rate  in  the  sample  and  f()  is  a  reference  value  of  this  rate  when  »0. 

When  i«  e  floes  hit  piers  under  natural  conditions,  the  strain  rate  decreases  with  the  distance 
from  the  pier  In  the  first  approximation,  only  the  plastic  zone  of  deformation  may  be  considered 
and  on  the  basis  of  ice  indentation  tests  performed  by  Korzhavin  this  zone  would  be  1.5  to  2.0  times 
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the  width  of  the  indenter,  here  BQ,  the  pier  width.  If  the  velocity  of  the  ice  floes  is  VQ  ft/sec  and 
the  length  of  the  plastic  zone  LQ  2 B0,  the  mean  strain  rate  in  that  zone  is: 


dL0/L0  VQ 


dt 


2  Br 


(17) 


For  the  same  pier  width  the  strain  rate  is  proportional  to  the  velocity  of  the  floe  in  the  range 
where  eq  15  would  be  applicable.  Extrapolating  to  higher  rates  of  stress  application,  Korzhavin 
proposed  the  final  formula: 


(18) 


where  VQ  is  a  reference  velocity  of  3.3  ft/sec  (1  m/sec). 

Because  other  experimenters  have  found  thi/  ice  strength  does  become  relatively  constant  for 
higher  rates  of  loading,’  11  17  a  range  not  tested  by  Korzhavin,  the  author  believes  that  this  correction 
on  the  strength  of  ice  should  not  be  taken  into  consideration  until  further  measurements  are  made. 

Influence  of  contact  between  the  flowing  ice  and  the  pier.  Because  of  the  brittle  fracture  of  ice 
floes,  there  is  never  a  complete  contact  between  the  pier  surface  and  the  ice.  The  effective  length 
of  the  contact  surface  between  the  ice  and  the  pier  can  be  taken  as: 

B  =  CB0  (19) 

where  B  is  the  effective  contact  width  and  £  is  a  reduction  coefficient  which,  according  to  Korzhavin, 
shall  be  taken  between  0.4  and  0.7  for  ice  of  medium  strength  as  in  Table  V. 

Table  V.  Values  of  contact  coefficient  £  according  to  Korzhavin. 

Pier  width  Velocity  of  ice  floes,  ft/sec 

(ft)  1.6  3.3  6.6 

10-17  0.70  0.60  0.50 

20-27  0.60  0.50  0.40 


Influence  of  the  form  of  the  indenter.  The  shape  of  the  cutting  edge  of  the  indenter  influences 
the  strength  of  ice  considerably  according  to  experimental  investigations  made  by  Korzhavin.  He 
studied  the  effect  of  six  different  indenter  shapes:  semicircular,  flat  and  wedge  shapes,  the  latter 
with  wedge  angles  2a  =  120°,  90°,  75°  and  60°.  The  tests  were  made  at  32°  F  (0°C)  and  the  rate 
of  indentation  was  1.5  cm/min.  Altogether  73  tests  were  made  of  which  52  were  with  lateral  confine¬ 
ment. 

The  experimental  results  can  be  summarized  as  follows: 

1.  The  shape  of  the  crushed  zone  near  the  indenter  depends  on  the  shape  of  the  indenter. 
For  flat  indenters  it  is  semicircular.  For  wedges,  the  size  of  the  zone  decreases  with  decreasing 
wedge  angles  and  it  disappears  completely  for  wedges  of  2a  =  60°. 
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The  pressure  on  the  ice  caused  by  tire  indenter  increases  continuously  but  irregularly 
with  peneirution.  At  the  instant  of  cracking,  the  pressure  falls  and  the  load  has  a  vibrating  charac¬ 
ter 

3.  The  indentation  of  ice  produces  mostly  a  brittle  failure  with  crack  formation. 

•1.  Tor  computing  the  pressure  of  large  ice  floes  on  piers,  Korzha'in  recommends  the  use 
of  a  shape  factor  whose  values  are  given  as  follows: 

a.  for  flat-edged  piers  (2u  180").  1.0. 

b.  for  semicircular-edged  piers,  £'  0.90 

e.  for  wedge-shaped  piers  (60°  ■  2«  120"),  C'  0.85  J  sin  a.  (20) 

Kor/huvm’s  formal  a  for  indentation.  The  basic  formula  (eq  11)  still  applies  if  the  crushing 
strength  of  ice  by  indentation  is  obtained  from  the  previous  considerations.  We  then  have: 
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2.5  <»0  for  B  15  B0. 


(21) 


When  the  kinetic  energy  of  the  ice  Hoe  is  not  high  enough  to  develop  the  maximum  pressure, 
relation  10  also  applies  but  in  that  case  it  must  be  recalled  that  o,  is  no  longer  a  constant  but  a 
function  of  the  pier  penetration  according  to  eq  21.  The  area  Uc  of  the  ice  sheet  required  to 
develop  the  maximum  pressure  can  be  determined.  For  a  triangular  nose  of  angle  2a,  the  integration 
of  eq  10  leads  to: 
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0.95  £  C  h  <>0 
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fan  a 


(22) 


winch  for  X()  B0  2  tan  a  d0  gives: 


im  V;. 


0  6  Hmaid0 


(23) 


where  <J()  is  still  the  length  of  the  cutting  edge  of  the  pier.  The  critical  size  of  the  ice  floe 
required  to  develop  this  maximum  pressure  is: 


(Mi  c  c‘  Bp,  • 


f>  Fp  tail  (i 


(24) 
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For  a  wide  floe.  B  >  1580,  this  is: 


1*5 C  C  ®o 
p  Vq  tan  a 


(25) 


Thus  for  average  values  of  the  parameters,  relatively  small  floes  of  100  >  100  ft  to  300  x 
300  ft  are  required  to  develop  the  maximum  ice  pressure.  The  force  that  a  smaller  floe  exerts  when 
stopping  in  front  of  the  pier  can  be  readily  found  from  eq  12.  For  a  triangular-nosed  pier,  it  is: 


Wmax  1.29  VQ  h  J  p  a,  tan  a  (26) 

which  is  not  much  different  from  Petrunichev's  formula  (eq  13). ’*  For  a  semicircular  pier  nose  the 
same  equation  can  be  used  in  which  Korzhavin  proposes  to  take  an  equivalent  wedge  of  angle  2a  = 
140°. 

Splitting  of  ice  floes  by  a  pier  with  a  vertical  edge 

We  will  now  consider  the  case  of  an  intermediate-size  floe  which  is  neither  stopped  by  the 
pier  nor  crushed  but  is  split  along  a  line  of  minimum  resistance.  The  minimum  resistance  may  be 
given  either  by  shear  at  an  angle  on  both  sides  of  the  pier  or  on  the  axis  of  the  pier,  or  by  tension 
cracks  forming  in  the  floe. 

Shear  cracks.  Let  us  consider  a  floe  of  size  L  x  fl  x  h  (Fig.  23)  whose  size  is  not  big  enough 
to  attain  the  critical  value  fic  given  in  eq  24.  The  shape  of  the  failure  crack  can  be  taken  as  a 
straight  line.  Contact  forces  between  the  pier  and  the  floe  can  be  reduced  to  two  normal  forces  T 
and  two  friction  forces  F.  For  equilibrium  conditions: 

P  -  2T  sin  a  -  2 F  cos  a  =  0.  (27) 

Because  we  have: 

F  I  T 

where  l  is  the  friction  coefficient  of  ice  on  the  pier  surface,  we  obtain: 

T  - 1 -  .  (28) 

2(sin  a  4  I  cos  a) 

Assuming  the  length  of  the  crack  Sj  we  get  along  the  line  of  shear  failure: 

T(cos  6  +  l  sin  0)  =  h  S  j  rQ  (29) 

where  rQ  is  the  average  splitting  stress  in  the  section  corresponding  to  the  shearing  strength  of  the 
ice.  It  is  evident  that  such  an  averaging  process  is  permissible  only  for  relatively  small  floes. 
Neglecting  friction  forces  in  the  first  approximation  we  get: 

2r0hS1sina 


sin  (a  4  0) 


(30) 


Figure  23.  Splitting  of  an  ice  floe  by  a  pier  having  a  verti¬ 
cal  edge. 40 


When  the  crack  reaches  the  side  of  the  floe,  we  have  S,  (B/2)  sin  fi.  The  force  P  becomes: 
rn(ifi  sin  a 

p  __r. — - .  (3i) 

sin  sin  (a  >  fi) 

The  splitting  occurs  for  Pmln.  when: 
fl  (90'  -  a). 

And  the  maximum  axial  force  acting  on  the  pier  with  this  lateral  splitting  is:  . 
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ilia  x 


2rn  B  h  tan  -  . 
0  2 


(32) 


If  the  floe  is  not  long  enough,  the  minimum  value  of  the  force  is  attained  with  a  longitudinal 
crack  in  the  axis  of  the  pier  and  this  force  is  then: 

H|nax  2r0L/isinu.  (33) 


Splitting  of  small  floes  by  tension  cracks.  For  small  ice  floes  and  very  sharp-edged  piers,  a 
splitting  failure  due  to  normal  tensile  stresses  at  the  crack  is  possible.  According  to  Korzhavin 
the  following  expression  is  suggested  for  that  case: 


/ i 
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max 


n  L  h  n  { 


(34) 


where  <»(  is  the  tensile  strength  of  ice  and  n  is  a  shape  factor  for  the  pier  varying  as  a  function  of 
the  wedge  angle  according  to  Table  VI, 
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Table  VI.  Shape  factor,  a,  aad  pier  angle,  2a.4* 

Pier  eagle  la  60°  70°  80°  90°  100 °  ISO0 

d  0.26  0.29  0.33  0.38  0.43  0.63 


In  all  cases  of  floe  splitting,  the  force  on  the  pier  is  always  smaller  than  the  maximum  force 
required  to  crush  the  ice  by  complete  indentation  through  the  floe.  Thus  in  the  case  of  large  floes 
the  splitting  is  always  partial  with  the  formation  of  cracks  of  limited  length. 

Fall  are  of  ice  floes  on  inclined  piers 

Inclining  the  cutting  edge  of  a  pier  produces  an  important  decrease  in  the  horizontal  component 
of  the  pressure  and  the  appearance  of  a  vertical  component.  This  facilitates  the  work  of  ice  cutting 
and  greatly  increases  the  stability  of  the  pier. 

After  a  floe  contacts  the  inclinSd  bdge  of  a  pier,  its  lower  sharp  edge  is  crushed  and  a  reaction 
appears,  as  shown  in  Figure  24,  which  can  be  decomposed  into  two  forces  V  and  T  as: 

V  =  T  cot/3.  (35) 

Under  the  action  of  these  components  the  floe  can  fail  in  three  different  ways: 

1.  by  bending  in  section  1-1 

2.  by  shearing  in  section  2-2 

3.  by  crushing  under  the  forces  T,  which  are  equal  to: 

H  -  2  T  sin  a .  (38) 

The  type  of  failure  that  takes  place  depends  oo  which  kind  of  stress  attains  its  ultimate  value 
first.  It  can  be  shown  thatjallure  by  shearing  is  the  most  probable  at  breakup.  As  is  known,  the 
ice  cover  of  rivers  and  lakes  represents  mainly  an  agglomerate  of  vertically  oriented  crystals, 
separated  from  one  another  by  narrow  layers  whore  impurities  can  collect.  During  the  spring  thaw 
the  resistance  of  ice  to  sliding  along  these  crystal  faces  is  much  reduced.  Observations  in  nature 
show  that  failure  by  shear  occurs  in  the  immediate  vicinity  of  the  pier,  so  that  a  channel  of  crushed 
ice  of  width  B0  remains  behind  the  pier. 

If  the  process  of  spring  thaw  is  not  yet  sufficiently  developed,  a  bending  failure  of  the  ice 
cover  is  possible.  In  that  case  too  the  fracture  occurs  in  the  vicinity  of  the  pier  at  a  distance  not 
exceeding  6  to  6  times  the  thickness  of  the  ice  cover.  The  width  of  channel  formed  by  the  pier  in 
the  ice  field  is  larger  than  in  the  previous  case,  but  is  still  only  a  little  larger  than  the  pier. 

Fracture  by  crushing  is  possible  only  for  rigid  piers  when  the  vertical  component  is  small. 

Failure  of  large  ice  lloes  by  shearing.  It  is  assumed  that  at  the  moment  of  shear  failure  the 
edge  of  the  floe  has  already  been  crushed  to  a  length  z  (Fig.  24).  Investigations  show  that  failure 
occurs  at  the  pier  surface  where  z  is  small.  We  then  find 

a)  The  length  of  the  shear  surface: 


(37) 
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Figure  24.  Failure  ol  ice  on  the  inclined 
edge  ol  a  pier.'0 
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and  the  calculated  forces  on  the  whole  pier  are  equal  to: 
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max 
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max 


In  the  case 


CB0hrQ 
sin  a 

1.1  C  B0  h  Tq  tan  /S 
sin  a 

of  a  semicircular  cutting  edge  we  get: 


^max  "  Y  <*B0h  ro 
Wmax  "  <BQhTQ  tan 


(42) 
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where  in  both  cases  the  contact  coefficient  C  can  be  taken  as  in  the  case  of  ice  crushing  (Table  V). 

Failure  of  large  ice  floes  by  bending.  This  case  might  theoretically  be  investigated  with  the 
theory  of  a  plate  lying  on  an  elastic  foundation,  as  the  inverse  problem  of  the  bearing  capacity  of 
an  ice  sheet.  This  theory  shows  that  the  failure  of  the  ice  sheet  begins  with  the  formation  of 
radial  cracks  starting  from  the  load  and  the  collapse  occurs  when  a  circumferential  crack  develops. 
This  latter  condition  gives  the  maximum  loading  condition  on  a  pier.  However,  observations  by 
Korzhavin  do  not  substantiate  the  theory  for  this  application. 

The  plate  theory  is  derived  for  vertical  loading  only.  Here  the  horizontal  component  of  the 
force  on  the  pier  combines  with  the  vertical  one  to  modify  the  stress  distribution.  The  theory  is 
applied  for  a  uniformly  distributed  load  over  a  circular  area.  Here  we  have  peripheral  loading 
along  different  shaped  pier  noses  and  observations  show  that  these  forms  are  a  prime  factor 
affecting  the  ice  pressure. 

In  elastic  deformation  the  theory  shows  that  the  circumferential  crack  would  appear  at  a  distance 
on  the  order  of  twenty  times  the  thickness  of  the  sheet  while  observations  by  Korzhavin  show  that 
the  dimensions  of  broken  pieces  of  ice  do  not  exceed  3  to  6  times  this  thickness,  and  are  frequently 
smaller.  In  the  theory  of  the  ideal  plastic  plate  the  position  of  the  circumferential  crack  has  to  be 
indirectly  assumed. 

We  will  thus  present  here  Korzhavin’s  semi-empirical  derivation  where  the  main  hypothesis  is 
that  the  collapse  of  the  ice  Bheet  occurs  at  a  distance  of  3h  from  the  sides  of  the  piers.  This 
distance  was  also  chosen  because  it  is  safe  for  engineering  design,  since  it  gives  a  higher  force 
than  could  probably  act. 

When  an  ice  floe  contacts  a  massive  pier,  indentation  occurs  first  when  the  lower  edge  of  the 
floe  is  crushed.  The  ice  then  comes  in  contact  with  the  side  of  the  pier,  a  vertical  and  horizontal 
reaction  is  set  up,  and  the  ice  is  ruptured  with  the  formation  of  at  least  three  radial  cracks  as 
shown  in  Figure  25.  These  radial  cracks  are  followed  by  transverse  cracks  breaking  a  piece  of 
the  floe  of  form  1-2-3-4.  After  that,  the  force  on  the  pier  is  reduced  until  a  new  contact  is  restored 
and  the  same  process  continues.  In  this  manner,  the  resulting  reaction  of  the  pier  on  the  ice  floe 
is  reduced  to  two  forces  on  each  side,  V  and  T  (Fig.  24). 


32 


ICE  PRESSURE  ON  ENGINEERING  STRUCTURES 


Figure  25.  Progressive  bending  tailures  of  an  ice  sheet 
at  an  inclined  edge  ol  a  pier. 40 


Again  taking  into  account  the  dynamic  friction  force  of  ice  against  the  structure  with  a 
coefficient  f  =  0.1,  we  get: 

H  =  2  T  sin  a  (46) 

H  .  2.2  V  tan  ft .  (47) 

As  the  separation  of  the  ice  into  two  strips  precedes  its  breaking,  the  greatest  pressure  on  the 
ice  occurs  at  the  moment  of  failure  along  the  line  3-4.  The  length  of  this  line  can  be  written: 


‘34 


"0 

2  sin  a 


(48) 


where  nQ  is  a  coefficient  higher  than  unity  taking  into  account  the  plan-shape  of  the  pier.  Korzhavln 
considers  the  strip  2-3-4-S  as  a  finite  element  clamped  in  section  3-4  and  acted  upon  by  forces  V  and 
T  at  the  other  end.  The  element  being  short  (3 h)  and  its  deformation  small  before  failure,  its  weight 
and  hydrostatic  reaction  are  neglected  in  relation  to  the  breaking  forces.  The  bending  moment  in 
section  3-4  is  then  given  by: 

U  =  V  %  -  T  y.  (49) 

Because  there  is  no  deep  indentation  of  the  ice  floe  (Fig.  24)  the  forces  V  and  T  are  not 
applied  at  the  center  of  the  ice  sheet  but  near  its  lower  surface.  Korrhavin  arbitrarily  proposed 
to  take  y  =  0.25  h  for  this  point  of  application.  At  failure  the  yield  moment  of  ice  for  an  elastic 
plate  is. 


(60) 


where  ob  is  the  flexural  strength  of  ice. 

With  these  relations  and  the  basic  hypotheses  x  =  3  A  and  y  =  0.25 h,  the  total  forces  acting 
on  the  pier  are  found  to  be: 
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=  0.9  Cq  h 
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max 


max 


(51) 

(52) 


and 


co  = 


0.73  n0 


12  sin  a  -  tan  /9 


(53) 


Assuming  a  similarity  between  the  bending  and  indentation  mechanisms,  Korzhavin  finds  a 
relation  between  the  coefficient  n0  and  the  angle  of  the  pier  2  a.  Adjusting  its  value  so  that  it 
will  fit  experimental  data,  he  finally  proposes: 

120° 

2.0 


2a  -  45° 

o 

o 

o 

75° 

90° 

105° 

nQ  -  0.94 

1.18 

1.42 

1.68 

1.98 

The  values  of  C0 

can  then  be  determined  from  eq  53  (Table  VII). 

Table  VII. 

Valiee  of  coefficient  C0. 

0 

2a  -  45° 

no° 

75° 

90° 

120° 

46° 

0.20 

0.17 

0.16 

0.16 

0.15 

60° 

0.24 

0.20 

0.19 

0. 18 

0.17 

O 

o 

r» 

0.38 

0.27 

0.23 

0.21 

0.19 

76° 

0.79 

0.38 

0.29 

0.26 

0.22 

It  can  be  seen  that  the  influence  of  the  angle  2  a  becomes  important  at  small  values  (2aC  60°). 
An  increase  of  Hmax  for  small  angles  2  a  can  be  explained  by  the  fact  that  a  sharper  pier  penetrates 
farther  into  the  floe  and  the  length  of  broken  pieces  is  correspondingly  increased.  If  the  cutting 
edge  is  rounded  in  plan,  some  additional  cracks  appear;  Korzhavin  then  recommends  taking  the 
rounded  edge  as  approximately  equivalent  to  a  sharp  edge  as  shown  in  Figure  26  whose  angle  a  is 
given  by  the  formula: 


2a  =  2at  +  —  (40°  -  at). 
B0 


(54) 


B©  2  “T* 


Figure  26.  Round  edge  and  equivalent  angle 
of  a  sharp-edged  pier.40 


i 
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Verification  and  utilization  of  method 

The  method  of  evaluation  of  the  force  of  impact  of  ice  floes  as  proposed  by  Korzhavin  has 
been  verified  In  nature  by  a  kinematic  method.  Knowing  the  mass  of  an  ice  floe  and  its  change 
in  velocity  during  its  action  on  a  structure,  it  is  possible  with  basic  laws  of  mechanics  to  deter¬ 
mine  the  reaction  of  the  structure.  This  was  done  by  taking  into  account  the  boundary  shear  force 
of  the  flow  and  the  component  of  weight  of  the  ice  floe  as  given  by  the  slope  of  the  free  surface  of 
the  river.  The  results  of  these  field  measurements  compared  to  computed  values  are  shown  in 
Table  VIII. 


Table  VIII.  Measured  sad  calculated  Ice  thrust.40 


Year  ot  measurement 


Data 

Unit 

1933 

1954 

1955 

1959 

1956 

1960 

No.  of  ice  floe 

1 

2 

3-5 

6-9 

10-11 

12-15 

Area  of  ice  floe 

1000  mJ 

15.6 

2.5 

24.0 

0.7 

3.6 

3.6 

Size  and  shape 
of  pier 

Width  ( B0) 

m 

4.4 

4.4 

4.4 

4.0 

7.7 

3.4 

Roundness  (radius  r) 

m 

1.0 

1.0 

1.0 

2.0 

0.0 

0.6 

Sharpness  (angle  2u) 

deg 

46 

46 

46 

80 

115 

83 

Inclination  (angle  fi) 

deg 

57 

57 

57 

45 

90 

90 

Measured  force  (Hmax) 

tons 

45.8 

17.5 

46.4 

15 

63 

28 

Calculated  force 

(Korzhavin) 

tons 

72.5 

25.5 

54.5 

27 

108 

70 

Thickness  of  ice  (h) 

m 

0.98 

0.50 

0.90 

0.85 

0.60 

0.70 

Note:  In  1954  and  1960,  the  tee  strength  was  lower  than  usual. 


The  calculated  values  are  in  general  greater  than  the  measured  ones.  It  is  always  very 
difficult  to  make  adequate  verification  in  the  field  because  of  the  extremely  variable  strength  of 
ice  during  the  spring  breakup.  To  our  knowledge  the  maximum  force  acting  on  a  pier  has  rarely 
been  measured  by  a  direct  method  in  nature.  With  a  dynamometric  plate  of  small  size,  Korzhavin 
found  a  maximum  of  130  psi  for  during  a  spring  breakup  and  the  Alberta  Research  Council,  in 
similar  circumstances,  found  a  maximum  pressure  ox  =  100  psi.  This  is  definitely  compatible 
with  design  values  of  ice  strength  suggested  by  Korzhavin.  For  such  plates  the  contact  coefficient 
would  probably  be  0.6  and  the  computed  design  indentation  strength  would  be  =  113  psi  (cf.  the 
measured  value  of  100  psi). 

To  summarize  the  results  of  the  Korzhavin  method  it  might  be  considered  that  the  force  exerted 
on  a  pier  is  the  smallest  of  that  necessary  to  break  the  ice  either  by  crushing,  bending  or  shearing 
as  computed  by  the  following  equations: 

1)  Indentation  failure  for  all  types  of  piers: 


// 


max 


2.5  £  £'  <>0  Bq  h 


(55) 


with  £  given  in  Table  V  and  O  on  page  36. 
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2)  Shearing  failure  for  an  inclined  pier: 


H 


max 


1.1 


CB0h 


r0 


tan  jS 
sin  a  ' 


(43) 


3)  Bending  failure  for  an  inclined  pier: 

Hmax  =  C0  °b  h  B0  tan  Z3  (49) 

with  £0  given  in  Table  VII. 

In  those  formulas,  i3  the  maximum  horizontal  force  in  pounds;  oQ,  ob  and  rQ  are  the 
crushing,  bending  and  shearing  strengths  in  psi;  h  and  BQ  are  the  ice  thickness  and  pier  width  in 
inches  and  finally  C  C,  C0  are  dimensionless  coefficients. 

In  the  case  of  an  indentation  failure  it  is  necessary  to  check  (eq  25)  to  see  that  a  floe  of  the 
required  size,  having  the  required  energy  to  attain  complete  indentation,  can  exist.  This  size 
is  limited  by  river  characteristics,  spacing  of  piers,  etc. 

We  emphasize  the  fact  that  sloping  a  pier  nose  is  the  most  effective  way  to  reduce  the  horizon¬ 
tal  force  caused  by  ice  impact.  For  a  pier  nose  with  a  horizontal  angle  2a  =  90°,  an  inclination 
of  /3  =  45°  will  reduce  the  horizontal  force  as  much  as  tenfold  with  the  usual  values  of  the 
parameters.  Furthermore,  the  vertical  component  of  this  force  then  contributes  to  the  stability  of 
the  pier.  Finally,  the  possibility  of  jamming  by  bigger  floes  stopping  and  arching  between  the  piers 
is  then  much  reduced. 

The  influence  of  the  horizontal  angle  of  the  pier  nose  is  not  straightforward.  For  a  vertical 
edge  pier  nose,  the  sharper  it  is,  the  more  effective  it  is  for  cutting  ice.  The  reduction  of  the  a;  -le 
from  2  a  =  120°  to  2  a  =  60°  would  reduce  the  horizontal  force  on  the  pier  by  about  25%.  However, 
the  effect  is  in  the  opposite  direction  for  an  inclined  pier  nose  and  a  sharpening  of  the  angle 
increases  the  force  on  the  pier  because  longer  strips  of  ice  have  then  to  be  broken  on  both  sides 
by  bending. 

Miscellaneous  effects  of  moving  ice 

Abrasive  action  of  ice.  When  an  ice  pack  is  moving  either  during  breakup  or  during  wintertime, 
the  floes  may  come  in  close  contact  with  wharves,  piers,  retaining  walls  or  dams  and  may  cause 
abrasion  of  the  structures  at  the  water  level.  Deep  cuts  in  a  wall,  on  the  order  of  a  foot,  have 
brought  about  the  collapse  of  a  wharf.1*  Piles  have  been  cut  in  less  than  five  or  six  hours. “ 
Recently,  during  a  spring  breakup,  concrete  walls  were  locally  eroded  to  a  depth  of  about  1  to  3 


The  erosion  of  river  shores  by  moving  ice  is  well  known  and  may  be  much  more  severe  in 
northern  countries  than  scouring  action  by  water  alone.  Even  rock  banks  are  extensively  eroded  by 
ice.  One  of  the  main  drawbacks  of  shore  protection  with  boulders,  artificial  blocks,  sand  bags,  etc. 
is  their  poor  resistance  to  ice  abrasion.  This  action  is  particularly  severe  on  concave  river  banks 
and  an  interesting  solution  has  been  used  recently  in  Minneapolis  on  the  Mississippi  River.  A  log 
boom  was  attached  at  the  upstream  end  of  a  bend  and  allowed  to  swing  free;  it  stabilized  at  a 
certain  distance  from  the  bank  because  of  the  spiral  motion  of  the  flow.  Ice  floes  tending  to  hit 
the  bank  were  deflected  by  the  boom,  whose  inertia,  together  with  the  inertia  of  the  water  behind  it, 
easily  guided  the  ice  out,  avoiding  the  bank  (Figure  27). 

The  maximum  local  pressure  that  ice  can  exert  while  moving  along  a  longitudinal  structure 
corresponds  to  local  crushing  of  the  ice.  The  normal  component  to  the  structure  face  can  be 
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Figure  27.  Ice  boom  deflecting  ice.  Mississippi  River  at  Elk  River. 


computed  from  the  previous  indentation  formula,  with  no  increase  in  this  case  for  the  size  of  the 
indenter: 

l,n  £°0-  (B6) 

The  longitudinal  pressure  of  the  moving  ice  is  the  dynamic  friction  force: 

'  =  %  t  (57) 


with  l  0.1. 

These  are  in  effect  localized  point  loadings  and  should  not  be  applied  simultaneously  over  the 
whole  wall  surface.  Because  of  their  repetitive  character,  abrasioo  results  in  the  long  run.  That 
is  why  it  is  recommedded  that  the  surfaces  of  structures  subjected  to  severe  ice  conditions  of  that 
type  be  hardened  either  by  using  special  surfacing  concrete  or  by  adding  extra  protection  such  as 
granite  or  steel  revetments. 

In  some  cases  the  design  of  low  slope  lateral  works  may  be  advised.  An  ice  foot  will  then 
form  and  remain  in  place  at  breakup  time  to  take  up  the  impact  of  moving  ice. 

Ice  accumulations.  A  moving  ice  field  may,  under  the  action  of  water  currents  or  wind»;hit  the 
shore  of  a  lake  or  river  and  any  engineering  work  located  there.  It  will  then  induce  important  ice 
accumulations  while  slowing  down,  and  dissipating  its  energy  by  breaking  and  accumulating  smaller 
ice  pieces  on  its  edge.  Figure  28  shows  an  accumulation  of  more  than  25  ft  on  the  shore  of  Lake 
Erie.  Accumulations  of  this  size  have  repeatedly  been  reported  in  the  literature  and  one  close  to 
100  ft  high”  has  been  observed 

These  accumulations  may  have  harmful  effects  on  installations  and  works  merely  because  of 
their  bulk  but  the  total  force  they  exert  is  usually  considerably  smaller  than  the  initial  impact  of  the 
ice  on  the  structure.  When  at  rest  the  force  they  exert  may  be  estimated  by  simulating  the  accumula¬ 
tion  of  a  granular  soil.  The  horizontal  pressure  at  any  point  above  the  water  surface  is  then  given 

bv  (Fig.  29V 
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Figure  28.  Accumulation  ol  ice  on  the  shore  ol  Lake  Erie 


Figure  29.  Horizontal  pressure  ol  ice  accumulation  against 
a  vertical  wall. 


p  =  K*y  y'  (i  -  e),  y  <  y0 


(58) 


where 

Ka  tan2  (45° 


,  the  coefficient  of  active  pressure  of  the  ice  pieces 


y  =  height  of  ice  on  top  of  the  considered  point 
<t>  =  angle  of  internal  friction  of  the  ice  pieces 
y'  =  unit  weight  of  ice 

e  =  porosity  of  the  ice  accumulation  (ratiocff  voids  to  total  volume). 
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Under  the  water  surface,  the  pressure  decreases  and  is  given  by: 

p  Ka(y0  >'■  (l  -  e>  (y  -  Vo)  (y  -  y‘)  0  -  e>) 

for  y  >  y0 


(59) 


where  y  is  the  unit  weight  of  water. 

The  difficulty  lies  in  assuming  reasonable  values  for  e  and  0 .  The  densest  ice  pack  with 
snow  and  slush  accumulation  between  the  solid  ice  pieces  would  not  have  an  e  value  sensibly 
smaller  than  0.3.  The  value  of  <j>  may  be  estimated  by  noticing  that  most  of  these  ica  mounds  have 
an  angle  of  repose  close  to  30°. 

Ride-up  ol  ice  on  bridge  piers.  On  a  pier  with  a  sloping  face  the  broken  floes  may  ride  up  the 
face  to  a  certain  level  before  falling  off.  If  there  is  insufficient  clearance  between  the  water  and 
the  bridge  deck,  they  may  occasionally  hit  the  superstructure  and  cause  damage. 

Because  a  floe  of  any  possible  form  pushed  by  a  large  icefield  may  occasionally  hit  a  pier  we 
have  to  consider  the  longest  floe  that  can  ride  up  the  face  of  the  pier  without  breaking.  It  was 
found0  that  a  floe  having  the  form  of  a  thin  beam,  simply  hingad  and  pushed  by  a  large  field,  would 
be  the  more  resistant  to  a  ride-up.  With  most  of  the  floe  length  out  of  water,  the  critical  length  was 
computed  from  equilibrium  conditions: 


L 


it 


(’ 


25° 


where  Mb  is  the  resisting  moment  per  unit  width  for  the  ultimate  strength  of  ice 


%h2 

~6~ 


(60) 


and  ac  is  the  critical  angle  at  which  the  floe  is  most  apt  to  break.  For  values  of  ob  of  65  psi  at 
breakup  the  critical  length  is  given  by: 

Lc  =  15.8  (61) 

where  Lc  and  h  are  in  feet.  When  this  length  is  laid  along  the  inclined  face  of  the  pier  it  gi/es  the 
maximum  possible  rise  before  '?:*  floe  falls  away. 


VERTICAL  FORCE8  EXERTED  BY  ICE  ON  HYDRAULIC  STRUCTURES 


Introduction 

Among  the  forces  exerted  by  ice  on  hydraulic  structures,  one  of  the  least  known,  and  perhaps 
most  often  neglected  in  design  practice,  is  the  vertical  force  transmitted  to  a  structure  by  an  ice 
sheet  during  water  level  fluctuations. 

The  adhesive  strength  of  ice  on  various  construction  materials  is  high  and  considerable 
moments  and  forces  have  to  be  resisted  by  the  walls  of  structures  following  a  variation  in  water 
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level.  With  a  rise  in  the  level,  piers,  caissons,  pile  groups,  well  casings  and  even  rectilinear  wharf 
sheetings  may  be  partly  lifted  off  their  foundations  and  damaged.  Observations  in  northern  Sweden 
have  shown  that  even  large  caissons  can  be  lifted  by  the  ice  and  one  was  raised  more  than  20  cm 
(8  in.).41  Damage  can  also  be  inflicted  on  revetments,  gates  or  protruding  parts  of  dams  and  intake 
works  frozen  in  the  ice  sheet,  following  the  usual  lowering  of  the  reservoir  level  during  winter 
operation. 

The  theoretical  determination  of  the  vertical  force  exerted  by  an  ice  sheet  is  a  difficult  task 
for  a  number  of  reasons, 

The  first  difficulty  arises  because  water  level  fluctuations  occur  frequently  during  the  winter. 
The  initial  failure  produced  in  a  thin  sheet  becomes  the  weak  point  where  subsequent  failures  will 
be  located.  This  time-dependent  phenomenon  defies  analysis  and  is  particularly  apparent  in  tidal 
waters  where  high  ice  caps  build  up  on  the  upright  surfaces  of  structures.  Fortunately  these  caps 
do  not  usually  have  extensive  horizontal  projections  and  the  computation  can  be  carried  out  as  with 
the  original  boundary  if  the  maximum  thickness  of  the  solid  ice  sheet  that  can  be  formed  with  a 
stable  water  level  can  be  predicted. 

Another  difficulty  is  our  limited  knowledge  of  the  basic  creep  properties  of  ice  of  various  types 
of  internal  structure  and  qualities  of  the  melt,  and  our  inability  to  predict  them  at  each  particular 
site. 

Finally,  not  one  of  the  smallest  problems  is  the  complexity  of  the  analysis  of  the  basv  be¬ 
havior  of  a  floating  ice  sheet  subject  to  very  slowly  applied  loads,  if  the  viscoelastic  properties 
of  the  ice  have  to  be  taken  into  account.  Many  approximate  assumptions  have  to  be  made  to  get 
usable  answers.  One  is  to  use  the  theory  of  elasticity  with  brittle  failure  of  ice.  The  other 
method  is  to  determine  the  maximum  loading  conditions  for  the  ultimate  strength  of  ice  on  the 
basis  of  plastic  theory.  In  this  latter  connection  it  might  be  said  that  analyses  pertaining  to 
the  similar  problem  of  the  bearing  capacity  of  an  ice  sheet  are  most  useful  and  have  been  applied 
to  this  case. 

Lifting  force  on  isolated  structures 

Elastic  behavior  under  concentrated  loading.  The  case  of  a  rising  ice  sheet  attached  to  a  pier 
or  a  pile  can  be  ideally  represented  by  the  similar  one  of  a  concentrated  load  applied  to  a  large 
floating  sheet  of  uniform  thickness. 

The  first  approach  to  the  solution  of  this  problem  is  the  use  of  the  theory  of  elastic  deforma¬ 
tion  of  an  infinite  plate  on  an  elastic  foundation.  This  solution  was  first  obtained  by  Hertz*'  for 
the  case  of  a  concentrated  load. 

The  equation  of  the  deflection  w  of  such  a  plate  is: 


where  q  is  the  external  load  applied  to  the  sheet,  s  is  the  foundation  modulus,  and  D  is  the  flexural 
rigidity  of  the  plate: 


12(1  -  v2) 
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in  which: 

E  =  moduius  of  elasticity  of  the  plate 
h  =  plate  thickness 

v  =  Poisson’s  ratio,  which  is  usually  taken  as  %  for  ice.” 
The  complete  solution  of  the  Hertz  problem  has  the  form: 


w 


PL2  .  r 
_ —  kei  - 
2nD  L 


where  L  is  the  characteristic  length  or  radius  of  relative  stiffness: 


(63) 


and  the  function  kei  is  a  Kelvin  function  of  zero  order.” 

At  the  origin  we  have  kei  (r/L)  =  -rr/4  and  the  deflection  under  the  load  becomes: 


w 

max 


P 

8  s  L2 


(64) 


The  distribution  of  the  unit  bending  moments  due  to  the  concentrated  load  is  given  by: 


The  radial  moments  become  negative  at  some  distance  from  the  load,  their  numerically  largest 
negative  value  being  about  -0.02 P  at  a  distance  r  close  to  1.9  L.  The  tangential  moments  are 
always  positive,  being  infinitely  large  at  the  origin  and  getting  asymptotically  equal  to  zero  at  an 
infinite  distance  from  the  center. 

If  the  concentrated  load  is  applied  ever  a  small  circular  area  of  radius  a,  the  moments  take  a 
finite  value  under  the  load.71  As  the  load  increases,  the  bending  stresses  below  it  become  equal 
to  the  flexural  strength  of  the  ice  and  the  sheet  begins  to  yield,  leading  to  radial  tension  cracks 
at  the  bottom  (Fig.  30). 

Circumferential  collapse  of  the  ice  sheet.  After  radial  tension  cracks  have  appeared  in  an  ice 
sheet  submitted  to  concentrated  loading  the  deflection  of  the  sheet  increases  quickly  until  the 
bending  stresses  along  a  circumferential  section  become  equal  to  the  flexural  strength  of  the  ice 
and  a  circumferential  tension  crack  is  formed  on  the  top  (Fig.  30).  The  ultimate  bearing  capacity 
or  collapse  load  of  the  sheet  is  then  reached.  One  first  method  of  analysis  is  to  consider  an 
elastic,  homogeneous,  isotropic  material  for  the  floating  ice  sheet  that  has  been  cracked  in  several 
radial  sections.  The  stress  distribution  in  such  a  sheet  has  been  completely  solved.50  The  load 
leading  to  the  circumferential  collapse  has  been  obtained  for  various  values  of  the  radius  a  of 
loading  and  a  Poisson’s  ratio  of  !,  for  ice  (Fig.  31). 
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Another  way  to  determine  the  theoretical  collapse  load  is  to  assume4*  that  the  ice  sheet  is  a 
thin,  rigid,  ideally  plastic  infinite  sheet  which  can,  without  cracking,  resist  a  full  plastic  limit 
bending  moment  per  unit  length.  <»0  ft^/4,  where  aQ  is  the  yield  strength  of  the  material  in 

tension.  At  very  low  rates  and  long  duration  of  stressing,  the  strain  capacity  of  ice  is  large,  as 
shown  by  creep  tests.  Meyerhof44  says,  therefore,  that  in  this  case  the  moment-rotation  capacity 
of  ice  in  bending  is  adequate  to  enable  the  limit  moment  MQ  (using  the  plastic  section  modulus  and 
creep  strength)  to  be  mobilized  at  the  critical  sections.  Hence,  the  collapse  load  of  floating  ice 
sheets  under  long-time  loading  can  probably  be  estimated  as  proposed  for  ideally  plastic  plates. 

The  method  of  analysis44  assumes  that  the  material  obeys  the  Tresca  yield  condition  and  the 
associated  flow  rule  where  the  shearing  strength  equals  one  half  the  yield  strength  in  tension. 

For  a  large  floating  ice  sheet,  cracked  radially  into  numerous  long  wedges,  the  total  collapse 
load  at  the  tip  is  found  to  be: 
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Equation  6ti  is  shown  in  Figure  31  for  various  values  of  a/L.  This  figure  can  he  used  either 
for  ductile  failure  with  the  limit  plastic  moment  MQ  under  long  term  loading  conditions  or  for 
brittle  failure  with  the  elastic  resisting  moment  Mb. 

Some  approaches  have  been  made14  to  take  into  account  the  time-dependent  deformation  of  an 
ice  sheet  by  making  hypotheses  as  to  the  rheological  behavior  of  ice.  Further  information  is  re¬ 
quired  however  on  the  creep  properties  of  various  ice  types  before  the  results  can  be  applied. 

Application  ol  analysis.  Figure  31  summarizes  the  results  of  computations  which  are  pertinent 
in  obtaining  the  maximum  vertical  force  that  ice  can  exert  on  an  isolated  construction.  The  upper 
curve  was  obtained  by  overestimating  the  radius  of  the  circle  of  collapse  and  is  on  the  safe  side  so 
that  ttie  maximum  force  can  safely  be  predicted  for  the  usual  range  of  small  structures.  The  main 
difficulty  in  applying  these  results  is,  here  again,  in  choosing  an  acceptable  value  for  the  flexural 
strength  of  ice.  Surveys  have  been  made  on  the  similar  question  of  the  collapse  loads  on  floating 
freshwater  ice  sheets  used  as  roads,  bridges  and  landings  by  pulp  and  paper  companies  in  Canada.10 
The  collapse  loads  varied  mainly  from  about  P  100  h*  to  P  =  500  ft  ,  where  P  is  in  pounds  and 
h  in  inches.  If  this  is  compared  with  the  plastic  theory  for  the  highest  values  and  the  elastic  one 
for  the  lowest  we  find  that  the  flexural  strength  of  ice  varies  from  ab  =  80  to  200  psi.  This  is 
compatible  with  results  obtained  directly  from  standard  flexural  tests  of  river  and  lake  ice."  J' 

The  case  which  has  been  considered  up  to  now  is  that  corresponding  to  a  load  distributed  over 
the  area  of  a  small  circle.  The  case  of  a  load  distributed  over  the  area  of  a  rectangle  may  be  ob¬ 
tained  by  methods  developed  for  simply-supported  plates.74  The  equivalent  of  a  square  area,  in 
particular,  is  a  circle  of  radius  r  -  0.57  a,  a  being  the  length  of  the  side  of  the  square.  The  effect 
of  any  group  of  concentrated  loads  on  the  deflection  of  the  infinitely  large  plate  can  be  calculated 
somewhat  by  summing  up  the  deflections  produced  by  each  load  separately.44 

Actually,  the  reactions  are  not  distributed  over  the  area  of  the  structure  but  along  its  periphery. 
It  is  thus  difficult  to  derive  the  stress  distribution  theoretically.  For  construction  of  a  rectangular 
form  it  may  be  safe  to  replace  the  shape  by  a  form  with  straight  sides  and  semicircular  ends.  Then 
the  lifting  force  on  the  straight  sides  can  be  computed  in  the  same  manner  as  for  a  long  straight 
wall. 

For  a  group  of  piles  of  rectangular  overall  form,  the  U.S.S.R.  standard  recommends  that  the 
total  uplift  force  be  taken  as  for  a  circle  of  radius  r  =  y  a b,  where  a  and  b  are  the  sides  of  the 
rectangle. 
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Verticil  force  on  a  straight  wall 

This  case  can  be  simulated  by  a  strip  load  acting  on  an  infinite  plate  on  an  elastic  foundation. 
If  a  number  of  equal,  central  concentrated  loads  are  applied  in  line  on  a  floating  ice  sheet,  the 
analysis  is  simplified  by  ignoring  longitudinal  bending  andu3ing  an  equivalent  strip  load  of  2  p 
per  unit  length  (Fig,  30). 

The  differential  equation  of  the  deflection  of  this  long  beam  is: 


diw  q  -  sw 
dx*  0 


(67) 


where,  as  previously,  q  is  the  applied  load,  s  the  foundation  modulus,  and  D  the  flexural  rigidity  of 
the  plate. 

In  this  case  of  loading,  2p  per  unit  length,  the  solution  is: 


For  the  elastic  state  the  maximum  positive  transverse  bending  moment  Wf,  per  unit  length, 
under  the  center  of  the  load  is: 


(69) 


The  maximum  negative  transverse  moment  occurs  at  a  distance  of  about  x  =  2.2  L  from  the  load 
and  is,  approximately,  =  -0.14pL.  The  maximum  deflection  is: 


and  the  width  of  the  deflection  band  is  about  x  =  3.3L.  The  yield  load  at  which  the  first  longi¬ 
tudinal  crack  develops  at  the  bottom  of  the  ice  sheet  below  the  load  can  be  estimated  from  eq  69. 
A  strip  load  has  to  be  more  than  3L  from  the  edge  of  large  sheets  to  satisfy  the  above  equations 
and  develop  the  estimated  yield  load. 

Chankin14  has  proposed  a  formula  very  like  eq  69  for  computing  the  maximum  force  acting  on 
a  straight  wall.  He  found  this  formula  to  be  valid  for  a  variation  of  water  level  of  from  0.3  to 
1  in. /hr,  ice  thickness  3  ft,  with  a  value  of  £  of  60  x  104  psi  as  obtained  from  measurements  in 
an  artificial  reservoir. 

In  the  case  of  a  very  rapid  rise  of  water  level  AH,  Lofquist42  proposed  to  use  the  same 
formula  when  the  ice  is  considered  to  be  an  elastic  material.  We  then  have: 


p  =  s  L  \H/2 


(71) 


where 

p  load  on  the  structure,  lb/ft 

s  =  foundation  modulus,  62.4  lb/ft3 
AH  water  rise,  ft 

L  -  radius  of  relative  stiffness,  ft. 
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This  nives  (he  unit  load  directly  as  a  function  of  the  variation  in  water  level.  In  contrast  with  the 
case  of  the  isolated  load  on  an  ice  sheet,  the  first  failure,  directly  in  contact  with  the  wall,  is  the 
final  one  as  the  wall  cannot  transmit  any  more  load  to  the  sheet  thereafter.  With  a  value  of  E  of 
60  •  lO4  psi,  the  deflection  wQ  at  failure  is  2.7  in.  for  a  1-ft-thick  sheet  of  200  psi  bending  strength, 
and  4.8  in.  for  a  3-ft -thick  sheet.  These  values  of  the  water  level  fluctuation  can  be  attained  quite 
rapidly  in  a  river. 

The  maximum  load  can  also  be  obtained  for  the  ideal  rigid  plastic  plate.  The  elastic  yield 
moment  of  ice  Mb  should,  however,  correspond  to  the  full  plastic  bending  moment  MQ.  The  only 
difference  is  the  change  of  Mb  by  Uu  and  the  maximum  possible  vertical  force  is  given  by  the 
same  equation  (eq  69)  which  gives  here  a  higher  value  for  p. 
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Strength  of  adhesion  of  Ice  to  structures 

In  some  cases  the  maximum  vertical  force  that  can  be  transmitted  to  a  structure  might  be  limited 
by  the  strength  of  adhesion  of  ice  to  the  material  the  structure  is  made  of.  Unfortunately  there  are 
very  few  available  data  on  this  basic  mechanical  property  of  ice  and,  as  for  other  properties,  there 
is  a  lot  of  scatter  in  the  experimental  results  because  of  insufficient  control  of  the  uniformity  of  the 
structure  of  the  ice  and  of  other  relevant  factors. 

It  has  been  shown17  that  the  adhesive  strength  of  ice  is  a  linear  function  of  the  temperature  up 
to  a  point  where  it  becomes  larger  than  the  internal  shear  or  cohesive  strength  of  the  ice  itself.  The 
sharp  transition  between  the  adhesive  and  the  cohesive  break  occurs  at  8V4°F  (-13°C)  for  strong 
snow  ice  between  two  stainless  steel  plates. 

Among  the  results  of  experimentation  for  artificial  ice  are  those  shown  in  Table  IX  which  were 
obtained  for  rather  high  rates  of  loading,  on  the  order  of  5  psi/sec,  and  results  depend  strongly  on 
the  rate.  It  must  be  noted  that  in  most  cases  the  strength  of  adhesion  of  ice  to  construction  materials 
is  as  high  as  the  shear  strength  of  the  river  or  lake  ice  itself,  which  is  usually  taken  to  be  from 
80  psi  to  150  psi. 7 
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Introduction 

Icebreakers  operate  either  under  arctic  conditions  or  in  ice  formed  in  a  temperate  climate. 

These  conditions  are  widely  different.  Arctic  ice  is  usually  several  years  old,  several  feet  thick, 
solid,  and  very  hard  because  of  the  low  prevailing  temperatures.  The  ice  in  the  temperate  zones 
is  softer  as  it  is  only  a  few  months  old  and  usually  not  more  than  3  ft  thick.  However,  it  has  the 
tendency  to  raft  and  form  pack  ice  of  extensive  dimensions.  In  rivers  there  are  often  accumulations 
of  frazil  and  slush  ice  right  down  to  the  bottom  of  the  navigation  channel.  We  shall  be  more  con¬ 
cerned  here  with  this  latter  type  of  ice  and  the  operation  of  icebreakers  in  navigable  rivers  and 
lakes. 

Icebreakers  are  best  defined  by  their  primary  functions;  to  break  solid  ice,  to  maneuver  in 
heavy  concentrations  cf  pack  ice,  and  to  clear  channels  through  which  other  ships  can  pass  in 
safety.  Ice-strengthened  ships,  in  contrast,  are  normally  cargo  transport  ships  strengthened  for 
use  in  ice.  Heavily  reinforced  hulls,  a  sloping  forefoot  and  a  ratio  of  horsepower  to  displacement 
of  more  than  unity  are  characteristics  common  to  both  types.  The  icebreaker  is  usually  distinguished 
by  more  powerful  engines,  higher  ratio  of  beam  to  length,  smaller  cargo  capacity  and  such  distinctive 
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Table  IX.  Adhesion  of  ice  to  various  materials  and  coatings  (Freiberger  and  Lacks”). 

Rate  of  load  application  is  5  psi/sec. 


Ice  adhesion,  shear  (psi) 

Metals: 

Mild  steel 

120 

Stainless  steel 

115 

Copper 

125 

Aluminum 

90 

Nickel 

85 

Zinc 

90 

Oalvanized  steel 

110 

Woods: 

Douglas  fir 

45 

Courbaril 

65 

Lignum  vitae 

80 

Plastics: 

Polymethyl  methacrylate 

40 

Cellulose  acetate 

25 

Glass: 

65 

Rubber: 

20  to  150 

Paints: 

Gray  deck  paint  (Navy) 

80 

Haze  gray  paint  (Navy) 

100 

Resin  films: 

Acrylics 

80  to  130 

Alkyds 

76  to  95 

Epoxy  s 

80  to  130 

Polyietrafluoroethylenes 

60  to  70 

Pol  y  tri  fluoro  chlo  to  e  Hi  y  1  en  e 

90 

Phenolics 

30  to  1 10 

Silicones 

40  to  95 

Urethanes 

95  to  130 

Vinylidenes 

50  to  80 

Vinyls 

90  to  11E 

features  as  heeling  tanks.  Icebreakers  designed  for  use  in  interior  waters,  as  in  the  Gulf  of  St. 
Lawrence,  may  possess  one  or  two  forward  propellers  which  are  not  used  in  polar  ioebreakers. 

The  first  ships  modified  for  breaking  ice  were  tugs  and  harbor  vessels  with  reinforced  bows  and 
increased  power.  They  cleared  ice  and  assisted  the  passage  of  cargp  ships  in  and  near  harbors.  The 
first  icebreaker,  in  the  true  sense,  was  Eisbrecher  1,  built  in  Hamburg  in  1871  for  service  on  the 
River  Elbe.  In  1890  the  Finnish  Murtaja  represented  the  greatest  advance  in  the  European  type  of 
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icebreaker.  It  was  able  to  break  20-in. -thick  solid  ice  with  no  snow  cover,  but  was  helpless  in 
ice-slush  and  snow-covered  ice.  The  full-lined  bow  pushed  the  snow  and  slush  forward  with  the 
result  that  the  ship  was  forced  to  ?  stop. 

At  the  same  time  the  development  of  icebreakers  on  the  Great  Lakes  proceeded  along  new 
lines  These  vessels,  which  got  into  difficulty  in  pack  ice,  were  able  to  force  their  way  through 
by  backing  into  the  ice.  The  natural  consequence  was  that  icebreakers  were  built  with  bow  propel¬ 
lers.  Thus  in  1888  and  1898  the  ferryboats  St.  Ignace  and  St.  Maria  were  built  with  both  stern  and 
bow  propellers.  The  primary  action  of  the  bow  propeller  was  to  wash  water  and  broken  ice  away 
from  the  forward  end  of  the  ship  and  thus  reduce  the  friction  between  the  ice  and  the  bow.  Further¬ 
more,  the  ship  could  break  into  pack  ice  by  alternately  running  backward  and  forward  on  the  bow 
propeller,  so  that  water  and  the  ship  were  forced  into  the  pack  ice.  The  work  of  one  bow  propeller 
is.  however,  asymmetrical  because  of  the  unbalanced  torque  reaction  and  in  modern  inland  ice¬ 
breakers  there  are  two  of  them,  rotating  in  opposite  directions. 

The  first  polar  icebreaker,  the  Russian  Ermak,  was  completed  in  Newcastle,  England,  in  1899. 
With  the  advance  of  shipbuilding  technology,  the  design,  power  and  size  of  icebreakers  developed 
steadily.  By  the  outbreak  of  World  War  II,  they  were  in  common  use  in  the  Baltic,  the  Gulf  of  St. 
Lawrence  and,  though  less  frequently,  in  the  Russian  and  Canadian  Arctic.  Military  demands  during 
the  war  stimulated  the  building  of  the  Russian  Stalin  class  and  the  United  States  Wind  class  ice¬ 
breakers.  Renewed  postwar  activity  has  included  the  building  of  increasingly  large  and  powerful 
icebreakers  such  as  the  Argentine  General  San  Martin,  the  Canadian  Labrador,  John  A.  MacDonald 
and  Louis  St.  Laurent,  the  Japanese  Fuji,  and  the  Russian  Moskva  class.  The  Russian  Lenin  with 
a  hitherto  unprecedented  nominal  horsepower  of  44,000  was  the  world's  first  nuclear-powered  sur¬ 
face  ship.  All  these  ships  are  primarily  polar  icebreakers.  The  development  of  other  icebreakers 
has  similarly  been  rapid,  particularly  for  use  in  the  Baltic. 

The  mechanics  of  icebreaking 

There  are  three  different  situations,  depending  upon  the  character  and  thickness  of  the  ice, 
which  an  icebreaker  may  have  to  deal  with:  the  breaking  of  a  continuous  ice  sheet  at  a  constant 
speed,  the  discontinuous  ramming  of  thick  strong  ice,  and  finally  the  movement  through  nonhomoge- 
neous  ice  consisting  of  a  mixture  of  rafted  ice  and  slush.  These  cases  will  be  dealt  with  succes¬ 
sively. 

Icebreaking  at  constant  speed.  An  icebreaker  acts  on  a  solid  ice  sheet  in  two  ways.  As  the 
vessel  moves  onto  the  ice  it  exerts  a  force  in  the  plane  of  the  ice  sheet  that  crushes  the  ice  in  the 
immediate  vicinity  of  the  stem  and  may  produce  long  radial  cracks  starting  from  this  point.  The 
forces  normal  to  the  ice  sheet  produce  bending  which,  initially,  causes  radial  cracks  to  emanate 
from  the  center  of  loading  because  of  tension  in  the  bottommost  fibers  of  the  ice  sheet.  As  loading 
is  increased,  the  radial  cracks  lengthen  and  increase  in  number  until  finally  a  circumferential  crack, 
roughly  at  right  angles  to  the  radial  cracks,  is  formed  by  tension  in  the  uppermost  fibers  of  the  ice 
sheet . 

The  icebreaking  operation  is  essentially  discontinuous  and  rhythmic.  Films  taken  during 
experiments  with  the  Mackinaw ”  and  other  enervations  show  that  the  ice  first  cracks  in  radial 
directions  from  the  stem  and  then  breaks  ahead  in  directions  at  right  angles  to  the  radial  cracks. 

The  broken  sectorial  pieces  turn  around  on  their  edges  and  pass  along  the  icebreaker’s  side.  Since 
the  overall  width  of  these  broken  floes  is  less  than  the  icebreaker’s  beam,  the  sides  of  the  bow 
come  in  contact  with  the  ice  sheet  and  a  further  two  or  three  rows  of  floes  are  produced  by  bending 
failure  Some  fleas  are  forced  under  the  ice  field  on  the  side  of  the  ship  bu'  the  majority  are 
broken  into  smaller  pieces  under  the  hull  and  appear  behind  the  ship  in  the  channel  it  has  made. 

This  breaking  action  is  essentially  of  the  vibratory  type.  When  the  bow  of  the  vessel  moves  forward 
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it  is  deflected  upward  over  unbroken  ice,  which  is  broken  up  on  the  downswing.  The  frequency  of 
the  vibrations  depends  on  the  length  of  the  broken  pieces  and  the  speed  of  the  ship.  The  amplitude 
builds  up  for  every  cycle  until  the  vibrations  get  out  of  phase  due  to  local  variations  in  the  thickness 
and  constitution  of  the  ice.  This  is  fortunate,  because  long-duration  resonance  vibration  could  easily 
lead  to  fatigue  failures  in  the  hull,  where  the  heaviest  bending  moment  occurs.  The  largest  ampli¬ 
tudes  measured  during  the  test  with  the  Mackinaw  were  8  mm  (Vlt  in.). 

For  the  same  power,  the  shape  of  the  bow  is  the  most  important  factor  affecting  the  efficiency 
of  the  icebreaking  operation.”  Because  the  flexural  strength  of  ice  is  normally  about  half  as  much 
as  its  compressive  strength,  the  general  concept,  of  course,  is  to  have  a  small  bow  angle  that  will 
increase  the  vertical  load  component  for  the  same  horizontal  one.  But  there  is  a  limit  to  the  sharp¬ 
ness  and  inclination  of  the  forefoot.  If  too  much  ice  is  broken  by  flexural  loading,  large  slabs  are 
forced  down  along  the  body  and  appear  at  the  propeller,  where  they  are  liable  to  impose  undesirable 
shock,  or  in  the  wake  where  they  interfere  with  successful  convoy  work.  If  the  bow  angle  is  steep 
and  the  buttock  lines  too  full,  the  resistance  set  up  by  the  broken  ice  as  it  passes  down  and  out  be¬ 
comes  great,  the  force  required  to  break  the  ice  in  compression  is  high,  and  the  efficiency  is  low 
for  the  same  power.  In  the  latter  stages  of  this  operation  the  ice  becomes  broken  so  small  that  it 
actually  looks  like  the  bow  wave  of  open  navigation.  In  the  final  instance  the  bow  should  be 
shaped  to  give  acceptable-size  floes  in  the  wake  with  the  minimum  expenditure  of  power. 

The  exact  shape  of  the  forebody  that  will  produce  such  a  result  in  different  cases  is  not  known 
because  the  force  that  each  type  needs  to  break  an  ice  sheet  of  given  thickness  and  quality  has  not 
been  determined.  This  is  so  because  of  the  extremely  complex  conditions  of  ice  loading,  the 
theoretical  difficulty  involved  in  the  analysis  of  the  stress  distribution  in  the  ice  sheet,  and  the 
little-known  behavior  of  the  ice  itself  at  collapse. 

A  physical  model  represents  the  best  approach  to  the  solution  of  that  problem  in  each  case  if 
a  suitable  material  can  be  found  to  represent  the  ice.  A  first  and  very  rough  approximation  can  be 
obtained  from  a  simplified  mathematical  model  where  only  the  vertical  component  of  the  force 
exerted  by  the  bow  is  taken  into  account.  The  loading  is  then  assumed  to  be  distributed  over  sectors 
of  a  semiinfinite  sheet  made  of  individual  wedges  resting  on  an  elastic  foundation  and,  finally,  the  . 
collapse  is  taken  to  be  that  of  a  perfectly  elastic  body.  The  analysis’0  then  shows  that  for  a  small 
area  of  loading  on  the  edge  of  a  semiinfinite  sheet  made  of  45°  unit  sectors,  the  vertical  force  at 
collapse  is  approximately: 

I 

P  =  3Mb  (73) 


where  *fb  is  the  maximum  yield  bending  moment  of  ice  per  unit  width 


o.  h 4 


(74) 


and  ob  is  the  flexural  strength  of  ice  and  /  its  thickness. 

Because  the  force  required  to  break  the  ice  sheet  depends  essentially  on  the  characteristics 
of  the  ice  and  the  form  of  the  bow,  this  most  important  resistance  that  ice  is  offering  to  an  ice¬ 
breaker  is  independent  of  its  speed. 

In  the  icebreaking  process,  it  is  recognized  that  there  are  two  types  of  ice  forces”  resisting 
a  ship  moving  in  a  solid  ice  sheet:  direct  ice  forces  and  dynamic  forces.  The  first,  which  do  not 
depend  on  the  icebreaker's  speed,  are:  resistance  to  breaking  of  the  ice  at  the  stem  and  along  the 
sides  of  the  ship,  resistance  to  tilting  and  submerging  of  the  floes,  dry  friction  of  the  floes  along 
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the  icebreaker's  hull,  ami  resistance  coming  from  wave  generation  during  the  oscillatory  moven^ept 
of  the  ship  and  the  floe.  The  dynamic  component  is  the  resistance  of  water  to  the  floe  movemapt 
to  winch  must  be  added  the  inertia  forces  which  occur  when  the  speed  of  the  pieces  of  ice  ip 
accelerated  from  zero  up  to  some  fraction  of  the  speed  of  the  vessel. 

The  direct  forces  depend  upon  the  characteristics  of  the  ice  andthe  dimensions  and  shapp  of 
the  icebreaker.  The  hydrodynamic  and  inertia  forces  also  depend  on  the  ice  dimensions  and  the 
form  of  the  ship  but  principally  they  depend  on  the  square  of  the  ship's  velocity,  if  the  spe^d  is 
not  too  low  The  following  seimempirical  formula*4  has  thus  been  proposed  to  represent  the  ice 
resistance  to  a  moving  ship: 

R  (C ,  .  C2  V 2)  Bh  .  (75) 

Where  R  is  the  ice  resistance  in  pounds,  V  is  the  speed  in  knots,  B  the  beam  in  f«et  and  h  the 
ice  thickness  in  inches. 

The  values  of  C  t  and  are  in  a  way  a  measure  of  the  efficiency  of  an  icebeakef  in  a  continuous 
operation  although  they  must  still  be  slightly  dependent  on  h.  There  are  still  very  fept  reliable  meas¬ 
urements  of  C(  and  C„.  They  are  of  the  order  of  40  and  1.5  for  the  Mackinaw 

Icebreaking  by  ramming.  For  very  thick  ice  the  icebreaking  process  is  no  longer  continuous  in 
the  forward  direction.  Instead  the  icebreaker  starts  at  some  distance  from  the  ice  and  runs  into  it 
at  maximum  speed.  As  the  speed  is  reduced  to  zero,  kinetic  energy  is  transformed  iqto  icebreaking 
work 

Because  all  available  energy  is  concentrated  on  breaking  the  thick  ice  sheet  by  bending  only, 
this  case  is  usually  considered  to  be  the  critical  one  in  the  design  of  polar  icebreakers.*'  In 
contrast  to  continuous  movement,  there  is  no  resistance  from  isolated  floes  in  the  channel,  and  in 
addition  to  the  engine  forces,  the  inertia  of  the  ship  is  available  to  break  the  ice 

The  fundamental  analysis  of  the  maximum  capability  of  a  ship  to  b^eak  ice  by  ramming  has 
been  developed  by  Vinogradov."  He  assumes  that  an  icebreaker  m  known  velocity  strikes 

a  uniform  ice  sheet.  The  bow  of  the  ship  glides  up  until  the  dov  ure  readies  a  magnitude 

which  causes  the  ice  sheet  to  collapse  as  can  be  estimated  rough.  q  73.  rle  assumes  that 

because  of  the  short  time  involved  the  thrust  is  horizontal  at  all  times  a  '  ‘  >at  t*  1  geometry  of  the 
water  plane  is  not  seriously  affected. 

The  principle  of  conservation  of  energy  is  applied.  The  energy  expended  is  the  ship's  kinetic 
energy  plus  the  work  done  by  the  propeller  thrust  acting  through  the  distance  traveled.  The  energy 
expended  is  directed  into  three  channels:  1)  energy  dissipated  by  impact  of  the  bow  of  the  ship; 

2)  potential  energy  due  to  the  ship  being  raised  and  changed  in  trim;  3)  frictional  logs  caused  by 
rubbing  of  the  ship  against  the  edge  of  the  ice.  In  the  case  when  the  ship  stops  at  the  instant  of  ice 
collapse 

E ,  .  E2  £.,  >  E4  *  E5  (76) 

where 

£,  kinetic  energy  of  the  ship  when  the  ice  is  first  touched 

£,,  =  energy  derived  from  propeller  thrust 

£3  energy  dissipated  by  impact 

£4  =  potential  energy  acquired  by  the  ship 

Eb  =  energy  loss  oy  friction. 
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L  Length  between  perpendiculars 
B  Beam  at  water  line 
D  Draft,  loaded 
V  Submerged  volume 
S  Draft  water  plane  area 
W  Displacement  =  loaded  weight 
a  Water  plane  area  coefficient  =  S/LB 
8  Block  coefficient  =  V/LBD 
G  Ship's  center  of  gravity 
CB  Center  of  buoyancy  center  of 
submerged  volume 


CF  Center  of  flotation  =  center  of 
water  plane  area 

Ml  Longitudinal  metacenter,  the  point 
about  which  the  ship  pitches 
mL  Longitudinal  metacentric  height 
GMl  (measures  stability) 

P  j  Maximum  value  of  vertical  reaction  at  stem 
T  Propeller  thrust 
4 >  Angle  of  stem  with  horizontal 
A  6  j  Maximum  angular  change  of  trim 
ADj  Maximum  change  in  draft  at  CF 


Figure  32.  Variables  and  nomenclature  used  in  ramming  analysis. 


Let  W  represent  the  weight  or  displacement  of  the  ship  and  K0  the  velocity  of  the  ship  when 
the  ice  is  first  touched.  Kinetic  energy  absorbed  during  the  operation  (assuming  negligible  speed 
after  impact)  is: 


S,  .  V*.  (77) 

The  next  item  considered  is  the  energy  delivered  by  the  propellers  to  the  ship  while  it  is 
gliding  up  on  the  ice  edge.  During  this  interval  there  is  a  reduction  in  mean  draft  designated  by 
ADj  and  the  ship  assumes  an  angle  of  trim  of  A0j  (Fig.  32).  The  distance  from  the  point  of 
contact  of  the  stem  to  the  center  of  flotation  is  designated  by  q.  The  stem  of  the  ship  is  sloped 
at  angle  4>  from  horizontal.  Then  from  the  instant  of  first  contact  until  the  ice  collapses,  the 
linear  advance  of  the  ship  is: 

A  D  j  cof  4>  +  Q  A  0  j  cot  4>  • 

Let  T  represent  the  average  value  of  propeller  thrust  during  this  advance;  then 

E2  =  T(AD  j  +  q  A0j)  cot  0  (78) 

This  formula  should  be  expressed  in  different  terms  so  as  to  include  P  v  the  maximum  value 
of  the  vertical  force  developed  at  the  stem.  Assuming  that  P1  is  small  compared  with  the  displace- 
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roent  W,  and  that  the  change  in  draft  and  trim  do  not  seriously  change  the  properties  of  the  water 
plane,  we  then  have  ADt  =  P  x/wx  S,  S  being  the  water  plane  area  and  wt  the  density  of  sea 
water.  The  angle  of  trim  A#t  depends  on  the  applied  moment  P^q  and  the  longitudinal  metacentric 
height  m:  thus  A0t  =  P^q/Wm.  This  energy  term  can  then  be  expressed  as: 


The  water  plane  area  S  equals  the  product  of  length,  beam  and  water  plane  area  coefficient, 

3  =  LB  a.  The  weight  of  the  ship  equals  the  product  of  length,  beam,  draft,  block  coefficient  and 
density  of  sea  water,  or  W  =  LBDSw New  nondimensional  coefficients  Jr  1  and  are  arbitrarily 
set  up  by  relationships  q  =  Jtj(L/2)  and  m  -  (Jfg  a2L2)/DS.  Substituting  these  new  quantities  in 
the  last  equation,  we  get: 

DP1 

£2  =  A— ±  T  cot  (79) 

where 


According  to  the  theory  of  impact,  when  two  bodies  collide  normally  there  is  always  a  dissipa¬ 
tion  of  energy,  whose  magnitude  depends  on  the  relative  velocity  ahd  a  physical  constant  e  known 
as  the  coefficient  of  restitution.  Now  the  stem  of  the  ship  does  not  collide  normally  withtfche  edge 
of  the  ice.  The  component  of  initial  velocity  VQ  which  is  directed  normal  to  the  edge  of  the  ice  is 
VQ  sin  \p  and  the  energy  dissipated  by  impact  should  be: 

Es  =  ^  <^0  8in<A)2(l  -  e2).  (80) 

The  author  does  not  believe  that  there  is  such  a  type  of  elastic  impact  as  given  by  Vinogradov 
and  this  term  should  instead  take  into  account  the  crushing  of  an  ice  wedge  in  the  ice  shelf 
corresponding  tn  the  horizontal  force  developed  during  the  prooess. 

/o,  y3  tan  ^ 

<7j  An  dn  =  -  tan  a  (81) 

u 

where  a,  is  the  indentation  strength  of  ice,  An  the  area  of  contact  between  hull  and  ice  and  dn  the 
distance  element  norma  to  the  shell  through  which  the  force  moves.  If  the  ice  deformation  is  small 
the  linear  advance  of  the  ship  corresponds  to  the  length  of  the  crushed  wedge  of  ice: 

y  =  ADjeot^  +  qAfljCot^ 

and  from  previous  considerations: 

A 3  D3P3 

Eg  =  Oq -  cot  i/f2  tan  a .  (82) 

3  W3 
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The  vertical  force  P  is  a  variable  which  keeps  increasing  as  the  ship  slides  up  on  the  ice. 

The  total  rise  of  the  point  on  the  stem  at  which  P  is  first  applied  equals  the  reduction  in  draft 
ADj  plus  the  angle  of  trim,  in  radians,  times  the  horizontal  arm  between  the  center  of  flotation 
and  the  stem  Adt  q.  The  potential  energy  set  up  by  the  force  P  is  therefore 

ADj  A0j 

E4  =  f  Pd  AD  +  f  P  q  d  Ad.  (83) 

0  0 

Energy  is  dissipated  by  sliding  friction  between  the  shell  plating  and  the  ice.  The  coefficient 
of  sliding  friction  /  must  be  applied  to  that  component  of  the  pressure  which  is  normal  to  the  plating. 
The  resultant  frictional  force  F  acts  in  a  direction  parallel  to  the  stem  of  the  ship  and  is  a  variable; 
half  of  it  acts  on  one  side  of  the  stem  and  half  on  the  other. 

The  energy  dissipated  by  frictional  force  F  acting  through  a  distance  determined  by  the  changes 
in  draft  and  trim  is  given  by: 


*5  = 


sin  if/ 


Ad, 

f 


Fd  AD  + 


sin  \l> 


a 

f  F  qd  Ad, 


(84) 


Consider  an  inclined  plane  intersecting  the  bow  of  the  ship  in  a  direction  normal  to  the  stem. 
This  section  of  the  bow  appears  as  a  wedge  with  essentially  flat  sides  and  the  normal  pressure 
ontthese  sides  makes  an  angle  /8  with  the  centerline  plane.  As  the  bow  rides  up  on  the  ice  shelf 
it  forms  a  wedgelike  groove.  Pressure  is  developed  normal  to  the  faces  of  the  groove,  and  friction 
aloag  the  faces  of  the  groove  directed  parallel  to  the  sloping  stem. 

Let  R  be  the  resultant  force  acting  normal  to  the  stem.  On  each  side,  then,  the  force  acting 
normal  to  the  plating  is: 


So  the  frictional  force  is  given  by 


The  magnitude  of  the  force  R  is  related  to  other  forces  acting  on  the  ship  as  follows: 
R  =  P  cob  \j/  +  T  sin  ifr. 


Equation  83  can  be  rewritten  as: 


F 


(86) 


Using  eq  84  in  eq  82,  combining  with  eq  81  and  substituting  the  previously  established  values, 
P  =  Wj  S  AD,  Pq  =  W  a  Ad,  and  others,  we  get: 


£4  ^ 


cot  *l>  \  **  1  + 
cos  /9/  2 


cob  /S 


AD 

T 


(87) 
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Substituting  all  the  values  of  component  energies  in  eq  83  there  finally  results: 

a  [7Pl  V  oy  Pi  T~ I  „  vo  [1  -  U  -  «2)  sin2 

LUy  W  W  J  ~  gD 


(88) 


in  which: 


1  - 

X  - 

r  f 
cos  B 

L  '  -i 

|  tan  y> 

1  + 

Y  - 

'  t  - 
cos  /9 

1 

cot  0 

1  + 

f  ‘ 
cos  /S 

cot  l/j 

cot  \p 


From  eq  76  we  obtain  the  downward  icebreaking  force  Pj  as  given  by  Vinogradov: 
Pl  =  XT  +  j*2T2  +  Y*-^-  [l  -  U  -  e2)  sin2  *  . 


(89) 


Taking  the  actual  energy  of  crushing  the  ice  instead  of  the  theoretical  energy  of  elastic  impact 
(eq  82)  we  get: 


s 

3 


where 


(90) 


S  = 


3W£ 


2  Oq  A2DzY  cot2  \p  tan  /S 


a  =  -  2XTS  -  — 
2 


ft  =  1  S3  +  2  *TS2  - 

27  3  ADg 


Tire  values  of  the  coefficients  X  and  Y  depend  on  the  dynamic  friction  of  ice  on  the  metal 
surface  which  is  usually  taken  as  0.10  to  0.15  for  fresh  and  brackish  water  ice  and  0.20  for  salt 
water  and  polar  ice."  One  of  the  important  factors  involved  is  the  presence  of  a  snow  cover  on  the 
ice  shelf  that  greatly  increases  the  coefficient  of  dynamic  friction  and  reduces  the  efficiency  of  the 
icebreaking  operation.  The  parameters  T,  V/,  VQ,  D  and  A  are  characteristics  of  the  ship  and  are 
discussed  at  length  by  Milano.47 
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Icebreaker  operations  in  pack  ice.  Accumulations  of  ice  floes,  snow,  slush,  brash  and  frazil 
ice,  usually  called  pack  ice  in  the  sea,  and  ice  dams  in  rivers,  often  produce  the  greatest  difficulties 
for  an  icebreaker.  It  is  not  possible  to  skirt  the  heaviest  positions  as  the  ice  conditions  underneath 
the  surface  are  usually  unpredictable. 

It  is  in  pack  ice  that  the  bow  propellers  of  icebreakers  are  exceptionally  useful.  They  produce 
water  lubrication  between  the  bow  of  the  vessel  and  the  ice.  Furthermore,  they  break  the  ice  into 
smaller  pieces  and  the  wake  of  the  icebreaker  becomes  clearer.  There  is  also  a  possibility  of 
working  on  the  ice  by  running  in  reverse  on  the  bow  machinery,  while  during  the  forward  motion  of 
the  stern  propeller  the  water  will  be  thrown  forward.  The  bow  propellers  should  be  inward  rotating 
so  that  the  wash  of  the  bow  is  more  powerful  and  effective  on  both  sides  of  the  vessel.  Observa¬ 
tions  have  shown  that  the  heaviest  mass  of  water  is  then  thrown  against  the  bow  at  approximately 
a  quarter  of  the  length  from  the  stem.  The  stern  propellers,  however,  should  be  outward  rotating  so 
as  to  throw  the  ice  pieces  outwards  and  give  a  clearer  track. 

Another  advantage  of  two  bow  propellers  is  the  better  maneuverability  they  give  the  vessel. 

At  low  speed  and  in  heavy  ice  the  action  of  the  rudder  is  weak,  but  by  suitable  regulation  of  revolu¬ 
tions  and  the  direction  of  rotation  of  the  propellers,  one  can  make  the  vessel  move  at  will,  even 
directlj  sideways.  But  the  bow-propeller  system  has  drawbacks.  The  costs  are  higher  than  with 
stern  p-opellers  only  and  the  water  and  ice  resistance  are  higher.  The  total  propulsion  efficiency 
is  also  smaller  for  the  same  power.  For  polar  operation  where  thick  hard  ice  is  encountered,  ice¬ 
breakers  are  not  designed  with  bow  propellers  because  of  the  decrease  in  ramming  ability. 

In  heavily  rafted  light  field-ice,  heavy  slush,  and  frazil  ice,  the  action  of  ramming  of  an  ice¬ 
breaker  is  not  really  efficient  because  the  ice  just  packs  closer  and  forms  a  stronger  cushion  after 
each  successive  charge.  One  practice  on  the  St.  Lawrence  River10  has  been  to  use  the  propellers 
as  drills  to  dislodge  the  ice  pieces  from  the  accumulation.  This  action  is  slow  but  positive  and 
sometimes  takes  place  under  conditions  where  the  ice  jam  is  so  high  that  parts  of  it  fall  on  the 
deck.  A  ship  with  two  forward  propellers  is  then  quite  an  asset.  Otherwise  the  ship  proceeds 
astern  with  its  stern  propellers  acting  as  drills.  All  Canadian  icebreakers6’  76  are  fitted  with  under¬ 
water  gear  which  will  stand  up  to  the  impact  of  breaking  ice  astern  and  propellers  which  will  in 
effect  drill  drift  ice  piled  as  high  as  20  to  30  ft. 

When  a  ship  is  moving  entirely  in  an  unconsolidated  ice  cover  or  pack  ice,  covering  part  or 
all  of  the  water  surface,  no  energy  is  required  to  break  the  solid  sheet,  but  the  friction  of  the  ice 
on  the  ship’s  hull  increases  considerably.  This  additional  friction  force  can  be  expressed  by  a 
formula  similar  to  the  one  used  to  determine  the  drag  resistance  of  a  ship  in  water  only:51 


R 


C  psv 
cd  - 


(91) 


where  R  is  the  ice  resistance  in  pounds,  V  the  ship’s  speed  in  feet  per  second,  S  the  wetted  surface 
of  the  hull  in  square  feet  and  p  the  density  of  water.  Figures  33  and  34  give  some  data  on  the  value 
of  R  for  icebreakers  of  the  U.S.S.R.  For  ice  concentrations  of  almost  100%  of  the  surface  area  an 
empirical  formula  has  been  proposed  by  Buzuev  and  Ryblin  based  on  measurements  on  the  Stalin 
and  Kapitan  Belousov:10 


R  =  8.6  L0,26  B0-25  (0.4  +  8.0  N pf26)  (92) 

where  R  is  the  ice  resistance  in  metric  tons,  L  and  B  the  ship  length  and  beam  in  meters  and  NFr 
the  Froude  number  of  the  ship  V /  /  gL  . 
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The  resistance  of  ice  pieces  accumulated  in  ice  ridges,  heavily  rafted  ice  or  frazil  and  slush 
ice  dams  is  unknown  at  present. 

Characteristics  of  icebreakers 

The  normal  task  of  an  icebreaker  is  to  assist  other  vessels  through  the  ice.  Other  tasks  include 
the  destruction  of  river  ice  jams  to  prevent  flooding  (Fig.  35),  winter  ferry  services,  supply  of  iso¬ 
lated  arctic  harbors,  patrol,  buoy  tending  and  hydrographic  surveying. 

When  assisting  vessels  in  smooth  waters  the  icebreaker  makes  an  open  track  and  the  other 
vessels  follow  as  closely  as  possible  without  risking  collisions.  At  sea  and  in  rivers  this  is  often 
difficult  because  the  track  does  not  stay  open  because  of  currents  and  wind.  One  or  more  of  the 
vessels  may  get  stuck  in  the  ice.  To  cut  it  loose  the  icebreaker  must  approach  closely.  Here  good 
maneuverability  is  of  great  importance.  The  ideal  is  to  have  two  icebreakers,  with  the  wider  of  the 
two  opening  the  way  while  the  other  offers  loosening  services  if  required,  as  done  on  the  St.  Law¬ 
rence  River.  In  very  difficult  conditions  it  may  even  be  necessary  for  the  icebreaker  to  tow  the  ship 
it  is  assisting 

The  characteristics  of  some  recent  icebreakers  are  shown  in  Table  X.  The  most  important 
main  dimension  of  the  icebreaker  is  its  beam,  which  determines  the  breadth  of  the  vessel  it  is  capa-. 
ble  of  assisting.  The  track  that  the  icebreaker  leaves  is  about  2  ft  wider  than  the  ship  itself.  With 
regard  to  length,  the  important  points  are  that  a  short  vessel  has  good  maneuverability,  while  a 
longer  one  has  more  room  for  machinery  and  engine  power  and  can  therefore  break  through  heavier 
ice.  For  the  majority  of  icebreakers  of  the  oceangoing  variety  the  usual  practice  has  been  to  adopt 
a  water  line  length/water  line  beam  ratio  in  the  range  of  4.0  to  4.5.  It  is  believed7*  that  there  is  an 
opportunity  for  improvement  in  polar  operations  by  increasing  the  length/beam  ratio  to  the  range  of 
5.0  to  5.5. 

If  there  are  no  restrictions  on  draft  due  to  water  depth,  the  draft  should  be  large  enough  to 
accommodate  large  propellers  placed  as  low  as  possible  so  as  not  to  draw  down  air.  As  the  stern 
propellers  are  usually  designed  for  higher  power  and  lower  revolutions  than  the  bow  propellers, 
their  diameter  is  larger.  This  difference  is  overcome  by  designing  the  icebreaker  for  a  2-ft  trim 
aft  which  also  makes  it  easier  to  reverse  out  of  ice  walls.  The  large  trim  and  tank  capacities  on 
icebreakers  make  it  possible,  within  certain  limits,  to  alter  the  draft  and  the  trim  of  the  vessel. 

The  typical  proportions  and  lines  of  an  icebreaker  (the  U.S.S.  Glacier )  are  shown  in  Figure  36. 
The  upper  portion  of  the  stem  is  nearly  vertical.  At  water  line  and  below,  the  stem  slopes  close  to 
30°  and  there  is  a  final  step  in  its  lower  portion  to  prevent  the  ship  from  riding  too  far  on  the  ice 
edge  and  getting  stuck.  The  midship  sectional  shape  is  close  to  a  circular  arc  with  a  center  above 
the  water  line.  The  virtue  of  the  circular  section  appears  to  be  that  it  facilitates  working  the  ship 
loose  by  means  of  the  heeling  tanks  when  beset.  The  first  effect  of  the  heeling  moment  is  to  break 
the  grip  of  the  static  friction  and  this  is  followed  by  a  loosening  of  the  ice  adjacent  to  the  hull. 

This  shape  tends  to  make  the  block  coefficient  very  small,  somewhat  lower  than  0.5  for  typical  ice¬ 
breakers.  It  produces  very  little  damping  to  any  induced  rolling  motion  and  icebreakers  are  heavy 
rollers  in  the  open  sea  if  stabilizing  fins  are  not  fitted.  The  after-body  form  of  an  icebreaker  usually 
narrows  to  a  point,  and  a  generous  clearance  is  provided  between  propeller  blades  and  shell  so  as 
not  to  let  any  ice  lodge  in  between. 

All  modern  icebreakers  have  heeling  tanks  that  can  be  put  to  use  in  one  or  two  minutes.  By 
pumping  water  from  one  side  to  the  other  a  heel  of  5  to  8°  to  either  side  can  be  achieved,  which  is 
normally  sufficient  to  loosen  the  vessel  from  the  ice.  The  capacity  of  ballast  tanks  in  icebreakers 
is  high  and  m  y  reach  15%  of  the  displacement.  In  certain  cases  the  icebreaking  efficiency  can  be 
increased  by  pumping  water  between  trim  tanks  in  the  fore  and  aft  peaks,  thus  changing  the  trim. 


Table  X.  Icebreakers  without  fore  propellers.  (Compiled  from  various  sources.) 
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Type  erf  machinery  Steam  Steam  Steam  Diesel  el.  Skinner  st.  Diesel  el.  Diesel  el.  Diesel  eL  Diesel  el.  Atomiosteam  turbine  Diesel  eL 


WL'Wottr  Lint  OWL »  Drof  t  Waft r  Lint  (B)' Dow  and  Buttock  Lint 

AP»Afttr  Ptrptndicular  FP»  Forward  Ptrptndiculor 


Main  Otck 


32' (WL) 


Draft  D:  25  ft  9  in. 
Displacement  at  DWL:  7400  tons 
Length  on  DWL,  L:  290  ft  0  in. 
Beam  on  DWL,  B:  72  ft  6  in. 
L/B  ratio:  4.00 


Block  coefficient:  0.48 
(Displaced  volume/LBD) 
Midship-section  coefficient:  0.77 

Freeboard  amidships:  12  ft  7  in. 


Figure  36.  Lines  and  particulars  ol  U.S.S.  Glacier. 
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Special  methods  are  sometimes  used  to  assist  in  the  icebreaking  process.  One  consists  of 
spraying  warm  water  from  the  circulating  pumps  onto  the  ice  at  the  stem.  Explosives,  salt  melting 
and  even  manual  ice  cutting  may  be  used  in  exceptional  circumstances.  Large  bow  vibrators  have 
been  used  in  small  river  icebreakers  in  the  U.S.S.R.  A  recent  ice  ferry  has  the  capability  of  ejecting 
warm  compressed  air  from  deeply  submerged  outlets  in  the  forebody. 

The  technique  of  icebreaking  in  Arctic  waters  consists  to  a  large  extent  of  searching  out  a  way 
through  leads  in  the  ice  where  resistancs  is  least.  Since  the  view  from  the  mast  is  comparatively 
limited  several  arctic  icebreakers  carry  one  or  two  helicopters  which  can  take  off  from  the  aft  deck 
for  reconnaissance. 


ICE  MODELS 

s'i 

Introduction 

In  the  last  few  years  knowledge  of  the  mechanical  behavior  of  ice  has  increased  to  the  point 
where  some  basic  laws  can  be  formulated  and  scaled  down  for  physical  models  to  study  engineering 
problems  caused  by  ice.  The  first  such  models  were  built  just  after  World  War  II  in  Canada,24  with 
paraffin  to  model  the  ice,  and  in  the  U.S.S.R.,44  with  weakened  ice  and  paraffin. 

The  technique  of  ice  modeling  has  evolved  along  two  completely  different  lines.  On  one  side 
studies  were  made  on  moving  ice  in  rivers  and  the  similitude  was  limited  to  representing  the  hydro- 
dynamic  conditions,  the  ice  itself  being  considered  a  rigid,  unbreakable  body.  On  the  other  side, 
models  were  built  to  measure  the  forces  caused  by  the  impact  of  ice  on  structures,  particularly  on 
icebreakers,  and  the  similitude  was  aimed  mainly  at  scaling  down  the  strength  properties  of  ice, 
somewhat  neglecting  other  dynamic  aspects. 

Laws  or  complete  similitude 

The  laws  of  similitude  of  an  ice  model  can  be  obtained  by  extending  the  standard  methods  of 
hydraulics.41  A  hydraulic  or  mechanical  model  is  said  to  be  in  complete  similitude  if  it  has  geometric, 
kinematic,  and  dynamic  similitude.  The  model  is  geometrically  similar  if  the  distances  between 
corresponding  points  of  the  model  and  the  prototype  are  in  the  same  ratio  for  all  points  in  the  model. 
This  can  be  expressed  as: 

t  A  |  .  (93) 

m 

Where  l  and  lm  are  the  distances  in  the  prototype  and  the  model  respectively  and  A  is  the 
geometric  scale. 

The  model  is  kinematically  similar  if,  at  corresponding  geometric  points  and  corresponding 
times,  velocities  and  accelerations  are  in  constant  ratios.  Thus 

v  v.vm  (94) 

a  a, a  (95) 

where  v.  v  are  prototype  and  model  velocities;  a,  am  are  prototype  and  model  accelerations; 
and  v,,  a,  are  the  velocity  and  acceleration  scale  ratios.  It  thus  comes  directly  from  previous  rela¬ 
tions  that  the  time  scale  t,  should  be  related  to  other  scales  by: 
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A  =  M*  (96) 

A  =  a, (2.  (97) 

Finally  a  model  is  dynamically  similar  if  at  all  geometrically  equivalent  points,  for  correspond¬ 
ing  times,  the  forces  in  the  model  and  the  prototype  appear  in  a  constant  ratio: 

F  =  F,Fm  (98) 

where  F,  Fm  are  the  prototype  and  model  forces  and  F,  is  the  force  scale.  This  last  condition 
is  sufficient  to  determine  all  other  scales  of  the  model.  If  all  the  forces  appear  in  the  same  ratio 
F*  for  corresponding  places  and  times,  it  is  clear  that  the  equations  of  mechanical  equilibrium  will 
be  the  same  everywhere  in  the  model  and  the  prototype.  All  terms  of  the  equations  for  the  model 
will  differ  by  the  product  of  the  same  constant  from  the  equation  of  the  prototype.  This  constant 
would  be  F*  for  the  equations  of  equilibrium  of  forces  and  A F,  for  work  and  energy  equations. 

Model  and  prototype  having  the  same  equations  will  thus  behave  in  exactly  the  same  dynamical 
manner. 

This  last  condition  leads  to  a  simple  method  for  obtaining  all  other  scale  values.  It  is  only 
necessary  to  express  each  of  the  forces  that  has  to  be  taken  into  account  in  the  ice  model  and  to 
keep  the  same  force  ratio  F,  for  each  of  them.  This  will  be  done  in  the  more  general  case  of 
complete  similitude  of  an  ice  model. 

Forces  ot  gravity  on  water,  ice  and  structures  in  water.  The  force  of  gravity  can  be  expressed 
by: 


F  =  p  g  I3  (99) 

where  p  is  the  density  of  the  materials  with  subscripts  1,  2,  3 . ;  g  is  the  acceleration  of 

gravity;  and  l  is  a  characteristic  length  of  the  material. 

As  it  is  not  feasible  to  change  g  on  a  model, 

F.  =  p.  A3  (100) 

P*  =  Pi*  =  Pz *  =  P3 . 

If  water  is  used  on  the  model  it  then  follows  that  all  other  substances  should  have  the  same 
density  that  they  have  in  the  prototype: 

p*  =  P 1*  =  Pz>  =  P3»  = 

So  that  eq  100  then  becomes: 

F,  =  A3  (101) 

which  is  the  basic  expression  of  the  force  scale. 

Forces  ot  inertia  ot  water,  ice  and  other  bodies  in  water.  The  inertial  force  corresponding  to 
any  change  in  velocity  in  water,  ice  and  other  bodies  as,  for  example,  waves,  current  changes, 
slowing  down  icebreakers,  etc.,  can  be  expressed  as: 
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dv 

dt 


where  dv/dt  is  the  total  acceleration  of  the  representative  mass.  This  gives: 


(102) 


F. 


A3 


v. 


f. 


and  with  eq  96: 


F,  = 


A2v2. 


Comparing  with  eq  101: 

v.  =  /T  (103) 

which  corresponds  to  the  Froude  similitude  where  the  Froude  number  NFr  =  V/Sfgi is  kept  the 
same  on  model  and  prototype.  It  can  be  easily  seen  that  the  Froude  similitude  comes  out  by 
comparing  only  the  force  of  inertia  with  the  force  of  gravity.  If  there  is  any  acceleration  whatsoever 
in  an  ice  model  it  has  to  be  scaled  down  with  the  Froude  similitude. 

Forces  of  water  friction.  The  water  friction  forces  on  different  boundaries,  including  ice,  can 
be  represented  by: 


(104) 


where  T  is  the  tangential  friction  force,  C  is  the  friction  coefficient  depending  on  relative  roughness 
t/l  and  Reynolds  number  NRe,  V  is  a  representative  velocity  outside  the  boundary  layer  and  I2  is 
the  surface  on  which  T  is  acting. 

For  rough  turbulent  flow,  as  is  most  generally  found  in  natural  rivers,  the  friction  coefficient  is 
independent  of  the  Reynolds  number.  To  get  the  same  coefficient  on  the  model  and  the  same  quad¬ 
ratic  law  it  is  necessary  to  have  complete  roughness  similarity  and  a  Reynolds  number  high  enough 
to  ensure  rough  turbulent  flow  at  the  reduced  scale.  For  river  models  there  is  thus  a  limitation  in 
the  choice  of  the  scale  that  should  be  big  enough  to  obtain  this  acceptable  Reynolds  number. 

When  a  body  is  towed  in  water  the  friction  force  also  has  the  same  form  as  eq  104  and  becomes 
independent  of  NRe  if  the  velocity  is  high  enough.  The  friction  force  in  the  developing  laminar 
boundary  layer  then  becomes  small  before  that  of  the  fully  developed  turbulent  layer  and  the  same 
reasoning  applies  as  for  the  previous  case,  In  all  cases  the  conditions  of  similitude  from  eq  104 
are. 


F  » 


A"”  V 


2 


which  gives  exactly  the  same  as  eq  103: 
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The  friction  force  is  also  in  complete  Froude  similitude  if  the  Reynolds  number  is  high  enough 
on  the  model. 

Forces  o(  water  pressure.  The  forces  of  water  pressure,  which  are  always  normal  to  a  surface, 
can  be  expressed  by 

N  p/2  (105) 

where  p  is  the  water  pressure  and  f 2  is  the  representative  surface.  In  similitude  this  gives: 

F,  p,  \2 

which,  with  eq  101,  gives  the  pressure  scale  ratio: 

P.  =  A.  (106) 

The  ratio  of  inertia  forces  to  pressure  forces  corresponds  to  the  well  known  Euler  number. 

Ice  friction  forces.  The  friction  of  ice  on  other  ice  or  on  foreign  substances  can  be  expressed 

by: 


T  =  (N 


(107) 


where  T  and  N  are  the  tangential  and  normal  forces  and  f  is  the  static  or  dynamic  friction  factor. 

It  can  readily  be  seen,  from  similitude  of  the  forces,  that  the  friction  factor  has  to  be  the  same  in 
the  model  as  in  the  prototype. 

Internal  ice  forces  at  collapse  of  an  ice  sheet.  If  brittle  failure  of  a  floating  ice  sheet  of 
thickness  h  is  considered  for  the  prototype  and  the  model  and  if  also  in  both  cases  the  failure  can 
be  predicted  from  elastic  theory,  three  types  of  forces  have  to  be  taken  into  account  corresponding 
to  possible  mode  j  of  failure: 

Shear  failure  forces. 

T  =  r0AfC,(®  £*  *)  (108) 

Crushing  failure  forces: 

N  -  oQ  h  l  j,  (109) 


Flexural  failure  forces: 


h2  C, 


(a  b  c  \ 
\  L’  L’  L’  / 


(110) 


where  T,  N  and  P  are  the  forces  producing  failure;  rQ,  t.Q  and  ob  are  the  shear,  compressive  and 
flexural  strength  of  ice;  f  ,  a,  b,  and  c  are  representative  lengths  (the  last  three  for  conditions  of 
loading);  v  is  Poisson’s  ratio;  and  L  is  a  characteristic  length  of  the  ice  sheet  that  is  related  to 
the  position  of  the  line  of  collapse: 

L  =  E  ^  —  (HD 

V  12Pg(l  -  v2) 
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Comparing  these  relations  with  the  basic  force  similitude  given  in  eq  101  we  then  get 


ro«  11  q» 


b* 


A 


(112) 


and 


L,  A  or  £.  A (1  -  (113) 

Scales  lor  complete  similitude.  The  complete  similitude  of  an  ice  model,  taking  into  account 
the  previous  forces,  is  thus  the  Froude  similitude.  Using  water  as  a  material  in  the  model,  the 
following  scales  have  to  be  taken: 

Geometric  scale  A  =  f  /{ 

m 


Velocity  scale 

v,  =  JT 

Acceleration  scale 

a*  =  1 

Time  scale 

t .  =  yr 

Force  scale 

F.  =  A3 

Mass  scale 

m,  =  A3 

Pressure  scale 

P.  =  A. 

The  model  ice  sheet  should 

then  have  the  following  characteristics: 

Thickness  scale 

h„  =  A 

Strength  scales 

r0*  =  A 

°Q*  =  A 

°b •  =  A 

Young's  modulus  scale 

II 

• 

CO 

Poisson’s  ratio 

Friction  factors 

v  1 

f  same  in  model  and  prototype, 
f  etc  1 

The  surface  characteristics  of  all  solid  material  represented  in  the  model  should  be  the  same, 
to  scale,  as  for  the  prototype. 

Other  negligible  forces  in  nature  are  surface  tension  forces.  Care  should  be  taken  that  they 
do  not  play  a  role  in  the  model  and  in  particular  the  model  ice  should  get  wet  as  the  real  ice  does. 
In  the  same  order  of  ideas,  the  model  should  be  big  enough  to  ensure  turbulent  flow  as  in  the 
prototype. 
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Practical  limitations  to  ic«  modeling 

The  first  difficulty  in  trying  to  make  an  ice  model  arises  from  the  extreme  variability  of  the 
properties  oi  '.he  natural  ice  to  be  represented,  Ice  moving  down  a  river  may  vary  in  form  from 
frazil  particles  a  few  millimeters  in  diameter,  agglomerating  together  to  form  very  loose  floes,  to 
solid  ice  pans  many  feet  in  length,  width  and  thickness  that  will  pack  closely  together  in  a  jam. 

At  the  same  place  in  a  river,  there  might  be  no  drifting  ice  during  one  winter  (all  ice  forming  and 
melting  in  place)  and  the  maximum  amount  the  upstream  stretch  can  ever  produce  the  following 
winter.  An  icebreaker  may  progress  in  a  solid  continuous  ice  sheet  but  most  often  it  will  go 
through  extremely  variable  conditions  of  hummocked  ice,  following  leads  going  through  brash  ice, 
ridges,  etc.  In  the  same  way,  the  strength  properties  of  ice  vary  within  wide  limits  even  for  the 
same  type  of  ice  and  even  from  hour  to  hour  for  the  same  piece  of  ice  if  temperature  and  radiation 
conditions  keep  changing. 

Thus  it  is  impossible  to  design  an  ice  model  that  would  reproduce  all  these  natural  conditions 
at  the  same  time.  Because  only  one  type  of  ice  will  normally  be  represented  by  the  model  some 
design  values  for  the  ice  properties  have  to  be  fixed  at  the  start  to  be  representative  of  the  con¬ 
ditions  to  be  studied.  At  times  extreme  values  will  be  taken,  as,  for  example,  in  choosing  high 
values  of  ice  strength  for  icebreaker  design.  At  other  times,  average  properties  can  be  chosen, 
as  for  the  size  of  the  model  ice  floes  to  represent  conditions  of  frazil,  slush  or  ice  jams  at  break¬ 
up. 

The  main  difficulty  in  trying  to  make  an  ice  model  is  to  find  a  suitable  material  to  model  ice 
at  a  reasonable  scale.  Up  to  now  four  materials  have  been  used  for  ice  modeling:  polyethylene, 
wood,  paraffin  and  ice.  Polyethylene  and  wood  are  used  in  small  pieces  to  represent  drifting  ice 
floes.  The  density,  roughness  and  average  form  of  the  drifting  ice  can  be  reproduced  but  it  is 
impossible  to  break  the  ice  pieces  in  the  model,  so  nene  of  its  internal  properties  are  modeled. 
Weakened  ice  is  produced  by  freezing  brine  very  quickly  with  a  method  developed  by  the  Arctic 
and  Antarctic  Institute  of  the  U.S.S.R.  The  fine-sized  ice  crystals  thus  produced  have  a  flexural 
strength  of  the  order  of  5  psi,  which  would  make  for  a  scale  of  '/„  if  compared  with  200-psi  strong 
ice  in  bending.  The  modulus  of  elasticity  and  the  shear  and  compressive  strengths  are  not  then 
to  the  same  scale,  which  makes  for  a  distortion  in  the  results  that  is  hard  to  appreciate.  Further¬ 
more,  it  is  a  considerable  undertaking  to  operate  a  big  model  under  cold  room  controlled  condi¬ 
tions.  Paraffin  and  model-making  wax  have  been  used  to  study  the  internal  behavior  of  an  ice 
sheet.  These  materials  are  very  strong,  with  a  flexural  strength  on  the  order  of  150  psi,  so  that 
theoretically  the  similitude  is  practically  impossible  at  an  acceptable  scale.  In  some  recent 
models,  additives  have  been  put  in  the  wax,  such  as  soap  powder,  coal  dust  and  kerosine,  con¬ 
siderably  reducing  the  strength  of  the  model  ice.  Such  methods  have  led  recently  to  suitable 
materials  to  model  ice  at  an  appropriate  scale.71 

Models  representing  only  the  hydrodynamic  behavior  of  ice 

In  hydrodynamic  models  the  conditions  of  internal  collapse  of  the  continuous  ice  sheet  are  not 
represented.  The  ice  is  essentially  made  of  individual  unbreakable  pieces  whose  hydrodynamic 
behavior  is  simulated.  The  Froude  similitude  applies  well  and  the  only  difficulty  in  choosing  model 
ice  is  to  represent  the  dimensions,  forms  and  surface  friction  effects  of  the  floes. 

One  of  the  very  first  studies  made  to  simulate  such  an  effect  was  carried  on  in  a  ship-towing 
tank  to  obtain  the  resistance  of  an  icebreaker  in  brash  ice  of  various  surface  concentrations. 17  The 
model  ice  was  paraffin-molded  in  hexagonal  pieces  and  tests  were  made  at  various  scales  to  deter¬ 
mine  the  coefficient  of  ice  resistance  as  shown  in  Figure  37.  It  can  be  seen  that  scales  ranging 
from  l4o  to  l/100  seem  to  be  adequate  for  that  purpose  if  the  Froude  number  is  higher  than  0.15. 
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Figure  37.  Coelficient  of  resistance  in  ice  (*/,0  covered)  lor  scale  models  ol  icebreak¬ 
er  Lenin. 


A  big  ice  model  was  built  to  represent  a  stretch  of  more  than  5  miles  of  the  St.  Lawrence  River 
at  scales  of  l/l40  horizontally  and  V„0  vertically.7’  The  aim  of  the  study  was  to  determine  the  effect 
on  ice  jams  in  Montreal  Harbour  if  new  land  were  reclaimed  for  Expo  67.  This  site  was  known  to  be 
subject  to  frazil  and  slush  ice  jams  up  to  about  30  ft  in  depth  with  water  level  rises  approximately 
20  to  25  ft  above  summer  levels.  The  frazil  floes,  slush  and  floes  were  represented  on  the  model  by 
small  polyethylene  grains  whose  surfaces  were  treated  to  get  a  wetting  effect  and  increase  the 
roughness.  The  model  was  built  according  to  standard  hydraulic  laboratory  procedures  for  fixed  bed 
models  and  the  bed  roughness  was  adjusted  for  clear  water  conditions.  The  model  was  then  fed  with 
ice  by  automatic  feeders  whose  rates  were  simulated  according  to  computed  values  in  nature.  The 
ice  runs  were  then  made  on  the  model,  simulating  natural  conditions,  and  the  results  of  water  level 
measurements  on  the  model  and  in  the  prototype  are  shown  in  Figure  38.  In  nature  the  jams  form 
during  cold  weather  and  there  is  always  some  bonding  by  thermal  exchanges  between  ice  floes  at 
the  water  surface.  On  the  model  there  is  no  cohesion  at  all  between  the  ice  pieces,  and  the  pack 
is  much  less  resistant  to  thrust  and  shoves  than  in  nature.  The  ice  accumulations  are  thus  bigger 
on  the  model,  the  head  losses  increase  and  the  water  levels  are  higher.  The  model  thus  errs  on  the 
safe  side  and  gives  the  worst  limiting  conditions. 

Another  example  is  a  model  built  to  study  the  flow  of  frazil  ice  on  a  river  in  Iceland."  An  intake 
to  a  power  plant  had  to  be  designed  that  would  get  water  free  of  ice  while  the  heavy  frazil  ice  run 
would  be  disposed  of  at  the  top  of  a  regulating  weir.  Polyethylene  foam  of  very  small  size  was  used 
with  success  to  represent  the  heavy  balls  of  aggregated  particles  of  frazil. 

Ice  models  are  most  useful  to  study  floods  caused  by  ice  during  the  spring  breakup  of  rivers. 

For  those  conditions  there  is  practically  no  thermal  influence  during  the  relatively  short  period  of 
jam  formation,  and  model  results  may  be  expected  to  give  a  good  approximation  of  the  truth. 

Figure  39  shows  such  a  model  in  operation  where  a  temporary  ice-delaying  structure  made  of  a  boom 
and  a  control  weir  had  to  be  studied.  Conditions  of  ice  accumulation  and  dry  jam  formation  were 
studied  in  the  model.  Model  ice  consisted  of  '4-in. -thick  by  l-in.-square  polyethylene  blocks,  surface 
treated,  which  simulated  the  average  characteristics  of  the  ice  floes  at  breakup  for  this  site  at  a 
scale  of  '/1J0 .  One  of  the  difficult  aspects  in  this  type  of  study  is  the  representation  of  the  solid 
ice  sheet  being  broken  up  and  pushed  back  under  the  th.ust  of  the  unconsolidated  ice  accumulation. 

On  this  model  the  problem  was  quite  local  and  the  solid  sheet  in  the  pond  formed  by  the  weir  was 
just  simulated  by  a  floating  wooden  boom.  This  boom  was  first  placed  at  the  upstream  end  of  the 
pond.  When  the  accumulation  became  unstable  and  was  shoved  under  the  boom  it  was  moved 
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Figure  38.  Water  level  and  discharge  lor  model  and  prototype. 


Figure  39.  Model  of  an  ice-retaining  boom,  Laval  University  (1967). 


artificially  down,  by  trial  and  error,  to  a  new  section  where  the  pack  ice  was  again  stable.  The 
final  retaining  place  was,  of  course,  the  position  of  the  designed  structures  at  the  top  of  the  weir. 
This  process  is  indeed  quite  crude  and  there  is  no  doubt  that  there  is  a  lot  of  room  for  improve¬ 
ment  ir.  the  techniques  of  ico  model  studies  of  this  type. 

Models  for  atudjdnf  the  collapse  of  an  ice  sheet 

Models  in  which  a  solid  ice  sheet  has  to  be  broken  have  been  used  up  to  now  practically  only 
for  the  design  of  icebreakers  and  bridge  piers.  The  technique  is  very  young  and  even  the  basic 
laws  of  similitude  are  being  questioned.71 

The  main  difficulty  here  is  to  use  a  material  that  would  simulate  all  properties  of  ice  at 
suitable  scales,  such  as  those  normally  used  in  model  ship  towing  tanks.  If  a  too-strong  model 
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material  is  ust*d,  such  as  untreated  paraffin  or  more  or  less  normal  ice,  a  distortion  has  to  be  intro¬ 
duced  to  various  scales  of  ice  properties  in  order  to  interpret  the  results.”  A  distortion  could 
also  be  introduced  only  in  the  scale  of  the  model  ice  thickness.  If  this  scale  is  given  bv 


hm  h.  h  (114) 

the  distortion  factor  D  is 

D  A-.  (115) 

h. 


It  can  readily  be  seen  from  eq  99,  102,  104,  and  105  that  the  forces  of  gravity  and  inertia  on 
the  ice  of  the  model  are  distorted  by  a  factor  D  but  that  the  friction  and  pressure  forces  on  ice  are 
correctly  scaled  down. 

For  certain  problems,  like  the  ramming  of  an  icebreaker  to  a  stop  on  a  solid  ice  sheet  where 
there  is  no  ice  movement  until  failure,  the  external  forces  of  gravity  and  inertia  of  the  ice  shelf  need 
not  be  taken  into  account  and  a  distortion  is  possible  in  the  thickness  scale  of  the  ice.  The  only 
conditions  of  similitude  left  are,  then,  the  similitude  of  crushing,  sliding  up  of  the  ship,  and  final 
bending  collapse  of  the  sheet.  Starting  with  the  most  important  one  (eq  110)  concerning  the  final 
collapse  we  get 

F,  „0,h‘t  (116) 

Choosing 

L.  A  (117) 

in  order  to  get  a  conservation  in  horizontal  length  for  collapse  of  the  sheet  so  that  all  geometrical 
properties  of  the  ship  stem  would  be  at  the  same  scale  in  plan  relative  to  failure  cracks. 

The  crushing  and  shearing  forces  are  also  reproduced  with  eq  108  and  109: 

F.  r0,Ah.  <7C.  A2  (118) 

so  that  the  final  scale  ratios  required  for  simulation  of  the  ice  sheet  in  that  case  would  be: 


Geometrical 

A 

Thickness 
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Strength 
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Because  in  the  total  energy  budget  for  this  process  the  energy  required  for  crushing  and 
shearing  the  ice  is  usually  very  small,  these  aspects  may  be  left  out  in  the  actual  simulation.  With 
the  model  material  properties  obm  and  £m,  it  is  then  possible  to  fix  the  scale  values  A  and  D  for 
each  application. 

For  other  phenomena  like  the  continuous  movement  of  an  icebreaker  in  an  ice  field,  or  the  impact 
of  a  floe  on  a  pier,  the  ice  pieces  are  quickly  accelerated  or  decelerated  and  the  inertia  of  the  ice 
cannot  be  neglected.  In  a  primary  analysis  no  particular  element  of  the  energy  budget  seems  negli¬ 
gible  and  complete  similitude  is  normally  required.  It  must,  however,  be  realized  that  no  model  ice 
will  reproduce  simultaneously  all  the  properties  of  real  ice.  The  geometry  of  secondary  failure  might 
be  quite  different  with  model  ice  because  of  the  usual  nonisotropic  behavior  of  real  ice  and  its 
weakness  in  shear  along  vertical  planes.  The  size  of  the  small  broken-up  pieces  in  a  model  may  thus 
be  distorted,  changing  the  overall  friction  effects,  as  for  an  icebreaker. 

Future  of  ice  modeling 

There  are  very  few  mechanical  problems  brought  up  by  ice  that  are  amenable  to  an  appropriate 
solution  by  mathematical  computations  only,  because  of  always-complex  boundary  conditions.  Such 
is  the  case  of  a  river  section  where  the  flow  of  water  and  ice  combine  to  fix  the  water  profile.  Such 
is  the  case  of  the  loading  characteristics  of  an  ice  sheet  rammed  by  an  icebreaker  or  of  a  bridge 
struck  by  a  moving  ice  jam.  In  most  of  these  cases,  physical  models  are  the  most  valuable  tools  for 
getting  acceptable  answers  to  the  questions  being  raised.  A  better  knowledge  of  basic  ice  properties 
combined  with  the  improvement  of  materials  used  to  simulate  ice  in  models  will  certainly  lead  to  a 
great  breakthrough  in  the  application  of  engineering  science  for  the  fight  against  river  ice  and  its 
use  for  beneficial  purposes. 
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